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C3 , Abstract 

We consider a third order non-autonomous ODE that arises as a model of fluid accu- 
mulation in a two dimensional thin-film flow driven by surface tension and gravity. With 
the appropriate matching conditions, the equation describes the inner structure of solutions 
X^/y \ around a stagnation point. In this paper we prove the existence of solutions that satisfy this 

/**s , problem. In order to prove the result we first transform the equation into a four dimensional 

"J ' dynamical system. In this setting the problem consists of finding heteroclinic connections 

that are the intersection of a two dimensional centre-stable manifold and a three-dimensional 
centre-unstable one. We then use a shooting argument that takes advantage of the information 
of the flow in the far-field, part of the analysis also requires the understanding of oscillatory 
solutions with large amplitude. The far-field is represented by invariant three-dimensional 
subspaces and the flow on them needs to be understood, most of the necessary results in this 
regard are obtained in [7|. This analysis focuses on the understanding of oscillatory solu- 
^ . tions and some results are used in the current proof, although the structure of oscillations is 

t ' somewhat more complicated. 
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1 Introduction 

In this paper is to prove the existence of solutions of 

f d 3 H 

that satisfy the following behaviour 



£ 2 + a )H 6 = 1, sel. (1.1) 



H^-L as |fl->oo. (1.2) 
l£l 3 

This equation has been deduced in [6] (see also [5]). It arises in a two dimensional model describing 
steady coating of a bumpy surface by a thin-film approximation. In particular (jl.ip results in the 
particular case that the motion of the fluid is driven by a balance of capillarity and gravity effects. 

In some regions the curvature of the substrate induces capillary forces of the same order of 
magnitude than the gravitational ones. In the steady regime the model describing such flows has 
the form (cf. (©)): 



^(A Q(s)+e ^r j =o ' £>o (L3) 

where we have neglected some non-relevant terms. The variable s stands for the arc-length that 
parametrises the substrate, and h is the height of the fluid over this surface. The parameter e is 
the ratio of the characteristic height of the fluid and the characteristic radius of curvature of the 
substrate. The function Q(s) describes the balance between gravitational and the capillary forces 
induced by the geometry of the substrate, it measures the tendency of the fluid to move in the 
tangential direction to the substrate as a result of the afore mentioned forces. 

If the function Q(s) has a constant sign, the motion of the fluid takes place always in the same 
direction. In such case, (|1.3|) can be approximated by the leading order term of <|1 .3)) : 



^(Q(s)h 3 )=0. (1.4) 

However, this approximation breaks down and it cannot be uniformly valid for arbitrary values 
of s if Q(s) changes sign. In such cases (|1.4[) predicts the onset of regions, where Q(s) is close 
to zero, with infinite height h, i.e. the fluid accumulates in those regions. As a consequence, the 
approximation (|1.4I) must be replaced by the model (|1.3[) . In the particular case in which in most 
of the substrate Q(s) is positive, but there exists a sufficiently small region (of size to be 
precise), where Q(s) = 0, a boundary layer analysis shows that, under suitable non-degeneracy 
conditions, the height of the fluid can be approximated by means of (jl.ip . the height of the fluid 
becoming of order e~ if . This asymptotic analysis shows also that the solutions of (jl.ip describing 
the stationary flows in those regions must satisfy (|1.2p . 

Equations similar to () 1 . 3[) where the main driving terms are the gravity and the curvature of 
the substrate, have been obtained, in a slightly different context, in [17] and [18] . This model 
can be obtained also as a particular case of the ones considered in [21] for specific choices of 
the parameters. See also |12] for a model that neglects gravity. Similar problems have been 
investigated in relation with industrial applications, such as the drainage of (metal) foams (e.g. 
[2"2"]). manufacture of lenses (e.g. [TJ] and JT3], although in the later case the effect of gravity can 
be neglected. These works offer numerical as well as formal (using perturbation methods) results. 

It is our aim to study the solutions of (|l.ip - ([1.2p rigorously. The main result of the paper is 
the following: 
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Theorem 1.1 For any nel there exists a solution of satisfying ifJ.ffj) . 



We sketch the main ideas in the proof of Theorem ll.il We first observe that the terms £ 2 H 3 
and 1 in give the natural scaling H ~ with $ ~ 1 (cf. I|1.2p ). which gives the leading 

order behaviour of (jl.ip . namely, 

lei -1 ^ + 1 + a ier a> ) # 3 = i , 



A change of variables with behaviour r ~ JL£|£|§ as |£| qq gives the dominant balance problem 
fd 3 <5> \ 

I —g- + 1 1 $ J ~ 1 as |^| — »■ +oo , $^las |r| -> oo . 
In such set of variables and r (|1.1[) becomes an autonomous dynamical system of the form 

^3- + 1 = ^ - (a - l)(cos0) 2 - F(9) , — = (cosO)^ , 

that we shall denote by (D), for the unknown ($, d$/dr, d 2 ^/dr 2 , 9) e R+ x M 2 x [— tt/2, tt/2]. 
Here 6* is defined by £ = tan and the function F is a linear combination of $ and its derivatives 
with coefficients that depend only on 9 and that vanish at 9 = ±tt/2. Thus, this system has the 
property that the three dimensional subspaccs {9 = ±7r/2} are invariant and the flow on them is 
described by the ODE 

d 3 $> 1 

T^ + i = W < L5 > 

The system associated to (|1.5[) for the unknown ($, d^/dr, d 2 $/dr 2 ) was studied in |7j. It has one 
single critical point, P s = (1,0,0) and therefore (D) has two critical points, p_ = (1, 0, 0, — 7r/2) 
and p+ = (1, 0, 0, 7r/2). Then, the solutions of (jl.ll) satisfying the matching conditions (|1.2[) 
correspond to solutions of (D) contained in the trajectories that connect the critical point p- as 
r — > — oo to p+ as t — > oo. Or equivalently, they are contained in heteroclinic orbits connecting 
these two critical points. 

The existence of a heteroclinic orbit for (D) is proved by means of a shooting argument in 
the direction of decreasing r. The shooting starts close to the invariant manifold {9 = tt/2} and 
the final argument will require information of the flow on the invariant manifold {9 = —tt/2}. 
For that reason we shall need the following information on p. 51) . First, that the critical point P s 
is hyperbolic and has a one-dimensional stable manifold and a two-dimensional stable manifold. 
Secondly, we proved in |7j that the only possible asymptotic behaviour of solutions on the stable 
manifold correspond to either 

lim $(t) = oo (1.6) 

or to 

lim $(t) = with r* > -oo . (1.7) 

t->(t„) + 



We shall also recall later that (|1.5p has a increasing Lyapunov function, and that this in particular 
guarantees the non-existence of periodic orbits. 

To start the shooting we first prove that there exists an invariant two-dimensional centre- 
stable manifold V+ locally defined near the point p + . All the trajectories associated to (D) whose 
starting initial data is contained in V+ converge to p + as r — > oo. We can parametrise the set 
of trajectories in V+ by means of one real parameter v taking values in some large interval. The 
behaviours (|1.6[) and (|1.7[) define two sets of values v. We prove that for very large values of v the 
corresponding trajectory satisfies (|1.6[) . On the contrary, if v is very negative we show that there 
exists a t* = tJ^v) such that (|1.7[) . 

It turns out that the sets of values v such that the corresponding trajectories satisfy either 
(|1.6[) or (|1.7[) are disjoint open sets. This implies, the existence of i/'s for which the corresponding 
trajectory does not satisfy neither ()1.6j) nor (|1.7|) . 
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The final step is to show that the trajectories associated to such v are globally defined inrel 
and that they satisfy 



d<I> 



dr 



dr 2 



<C, for any r € K (li 



for some C > 0, and that 



lim 0(t)=~. (1.9) 

r— > — oo ^ 

The idea is that if (|1.8[) and (jl.9p are satisfied we can use the fact that the dynamics of (D) 
become close to the ones associated to the trajectories contained in the unstable manifold of P s 
for the system associated to (|1.5|) and the trajectories have no alternative but to approach p_ as 

T — > — OO. 

The most technical part of the paper is the proofs of (|1.8[) and of ()1.9j) . These require to show 
that oscillatory behaviours with large amplitude for the solutions of (D) as r — > — oo must have 
a decreasing amplitude for decreasing r if neither (|1.6j) nor (|1 .T[) take place. The key point is 
that the structure of oscillatory solutions can be identified by looking at the several asymptotic 
regimes of (| 1 . 1 [) . There are, in particular, two very distinctive ones. For instance, the balance 
£ 2 / 3 H ~ oo for very negative £ will be relevant in our analysis. In this case the behaviour of 
solutions is described by 

This equation can be integrated giving that, in such regions, H behaves like a fifth order polyno- 
mial. The solutions of (|1.10|) are in fact a two-parameter family of polynomials, as we shall see. 
On the other hand, if £ 2 / 3 i7 ~ on bounded intervals, the dominant balance there is given by the 
equation 

d 3 H 1 

W = hs- (Ln) 

The analysis of (jl.lip plays a crucial role in our proofs and was already studied in UJ. The 
possibility of alternating regions where either (|1.10[) or p. lip dominates, builds up a scenario 
where solutions with large oscillations exist: The bouncing region of the oscillations are described 
by (jl.lip and the maximum amplitude regions are close to solutions of (jl.lOp . This phenomenon 
has been already observed for (|1.5p in [53] and explored rigorously in [7J. 

In order to prove (|1.6[) and (jl.7p we exploit this mechanism of oscillation. We argue by 
contradiction and assume first that (|1.8p does not hold. This gives (after a number of technical 
lemmas) that there exists a sequence {r*} with lim^oo r* = — oo such that 3>(t*) is a local 
maximum and linin^oo $(t*) — oo. We use that the oscillatory solutions with very large amplitude 
for very negative values of r can be approximated, after a suitable rescaling, by a sequence of 
functions \£\ 2/3 n n {0 where each U n solves (fTTTOl) in intervals [f(r£$f), ?(C m )]- The values r™ n 
being such that $(r™ m ) is the minimum in (r*,r*_ 1 ). In particular, in such intervals % n (£) 
are close to a fifth order polynomial solving ([l.lOp . The matching between two consecutive such 
functions is done into the inner region where <E> and Hn become close to 0, as it turns out, this 
inner regions lies around r™ m . As we have mentioned the dynamics in such bouncing region are 
dominated by (jl.lip and the rigorous matching can be adapted from that performed for p.5p (see 
[7j): the study of p. lip reduces to the one of a phase-plane analysis in which the bouncing can 
be encoded into the behaviour of a separatrix. This object attracts trajectories for increasing £, 
implying that its behaviour is generic. Reading off this behaviour into the functions T-L n implies 
that in the outer region they behave as a polynomial with a double zero near £(r™"). This 
in particular reduces the family of polynomials that give the outer region around each t* to a 
one-parameter family. Moreover, this analysis allows to get information on the relative size of 
consecutive maxima and minima, namely that the sequence of the maximum values decreases and 
that the sequence of minimum values increases (as n — > oo) and these contradict the assumption 
that (H]) does not hold. 
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The paper is organised as follows. Section [5] is divided in three preliminary parts. First in 
Section 12.11 we give some results concerning (|1.5I) , most of which are proved in [7J . In Section 12.21 
we reformulate (jl.ll) as a four dimensional dynamical system and reformulate Theorem 1 1.1 1 in this 
setting. The third part is Section [2 . 31 where we prove the existence of the centre-unstable manifold 
around p + . Section [3] is devoted to the analysis of the behaviours (jl.6[) and (|1.7|) for (D); in 
Section 13.11 we show stability under small perturbations of solutions that satisfy either of these 
properties, and in Section 13.21 we give necessary conditions on solutions of (|1.1[) to satisfy cither 
(|1.6|) or (|1.7|) . With the analysis carried out up to here we can then prove in Section |4] that there 
exist solutions on V+ that do not satisfy neither (|1.6p nor (|1 . T[) . We continue by proving that these 
trajectories of V+ do satisfy (( 1 . 8|) and (|1.9[) . In order to do that we first find in Section [5] that if 
(|1.8|) is not satisfied the sequences {t*} and {t™™}, described above, are well defined. Second, 
in Section [S] we find the contradictory results that {$(r*)} is decreasing and that {3>(T™ n )} is 
increasing. This part is very technical and needs by itself a few steps. Thus, in Section 16.11 we 
identify the scales of the outer region and the approximating polynomials near local maxima. 
This is based on the analysis of the solutions of (|1.10[) that is carried out in Appendix [Bj In 
Section 16.21 we perform the right scaling of the solutions under consideration and identify the 
range in which they are approximated by the polynomials. In this section we also prove that the 
approximating polynomials must have a double zero. This step requires the analysis of (jl.lljl 
given in Appendix [Al as well as the matching lemma given in Appendix [Cl fa result that has been 
adapted from [7J). In Section [6. 3[ with detailed information of the matching regions, we derive an 
(iterative) expression that relates the elements of the sequence of local maxima and another that 
relates the local minima, and that contradict that (|1.8I) is not satisfied. Finally, in Section [7J we 
finish the prove of Theorem 11.11 

Finally, we recall that equations similar to (|1.5I) have been studied intensively, see [3], [5], 
M' EL dUi [I7J, [IB], [H] and [13]), to mention a few, where similar equations arise in several 
related physical situations. Rigorous results concerning such equations can be found also in [I] 
and, concerning travelling wave solutions, in [2], [3], [M], [15] and [20]. It is interesting to note 
that many of these models yield higher order ODEs describing oscillatory fluid interfaces. We 
refer to [7J, where this aspect and related works are put into context. 



2 Preliminaries 

2.1 A summary of results for (11.5ft 

We now summarise some properties of (|1.5[) . most of which have been proved in [7j and will be 
used later in the proof of Theorem 11.11 It is convenient to rewrite (|1.5I) in the equivalent form 

d$ dW 1 , . 

dr dr dr <J> 3 

we then have the following result. 

Proposition 2.1 (i) There is a unique critical point for 12. 1\) in the domain {$ > , W € 
R , * 6 K} given by: 

P S = ($,W,V) = (1,0,0). 

(ii) The point P s is hyperbolic, the stable manifold of il.5\) at the point P s is tangent to the 
vector: 

3-1 

-3"3 

and the corresponding eigenvalue is Ai := —3^. 
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(Hi) At P s there is a two-dimensional unstable manifold locally spanned by the eigenvectors V2 := 
(— 3s/6, 33/6, 1) T and V3 := (3s/6, 3e/2, 0) T . The eigenvalues associated to the plane 
spanned by {1*2, ^3} are A2 := 3^(1 + i3?)/2 and A3 = A2. 

(iv) The trajectories associated to (2.1)) that are contained in the stable manifold and satisfy 
($, W, ^ P s , behave in one of the two following ways for decreasing r: Either they are 
defined for all r gl and satisfy 

lim ($,W,*) = (00,-00,00) (2.2) 

r— > — OO 

or ? alternatively, there exists a r* > — oo suc/i i/iai 

lim $(r)=0. (2.3) 



Moreover, the points of the stable manifold associated to P s with $ > 1 satisfy \2. 2)) and 
those with <f> < 1 satisfy (2.3)) . 

(v) Suppose that there exist tq G K and Co > 1 such that 

(*(t), W(t), tf (t)) G j (*, W, *) G M 3 : 77- < $ < Co, -Co < < C" , -C < * < C c 

L L-0 

/or aZZ r < r 0; i/ien 

lim (#(r), W(r), *(t)) = P, (2.4) 

and £/ie corresponding trajectory is contained in the unstable manifold of P s . 

Proof. All the statements of this proposition have been already proved in [7] except for (v). In 
order to prove this, we use an argument similar to the one used to prove Lemma 2.4 in [7]. We 
first recall that there exists an increasing Lyapunov functional E associated to ()1.5[) : 

E := +7^7 + $ , ^ = * 2 >0. 
2$ 2 dr 

This and the assumptions made imply that 

* 2 (s)ds< 00. (2.5) 

-00 

Using (|2.ip . it then follows that lim T _ 5 ._ 00 ^(r) = 0. Indeed, arguing by contradiction, one can 
construct a sequence r„ — > —00 such that there exits a £0 > such that either ^(r n ) > Eq or 
*(tVi) < -£o- Then (f2~Tj) implies that d^ /dr > -1, so cither *(t) > £ + (t„ - r) for r > r„ or 
*(t) < -£o+(t„-t) for t > r n . But this contradicts (|23)l ) since cither / r T " +eo/2 (*(T)) 2 dT > £ 3 /8 

° r f T :-s o/ 2^^ 2dT ^ e3 / 8 for a11 n - 

Now the second equation in (|2.1[) implies that remains approximately constant as r — > —00 
in any finite interval of arbitrary fixed length L. Therefore, if there is a subsequence {t„} with 
linin^oo r n = —00 satisfying lim„^oo W(r n ) ^ 0, we obtain that inf rg [ Tn . Tn+ £] | W(t)| > Eq > for 
n sufficiently large. It then follows from the first equation in (|2.1[) that the condition ^ < $ < Co 
fails if L is assumed to be sufficiently large (integration on the interval (r„, t„ + £) for sufficiently 
large n gives that |$(r„ +£) - 3>(t„)| > £ i > 0, but |$(r„ +L) - $(r„)| < (C5 2 - 1)/C for all 
n and L). Therefore lmv^-oo W(t) = 0. 

Using the last equation in (|2.ip as well as the fact that lim T _ ) ._ 00 ^(r) = it then follows in a 
similar way that lmv^-oc $(r) = 1. This gives (|2.4[) and the result follows. ■ 

The next lemma gives the detailed asymptotic behaviour in both cases (|2.2p and (|2.3j) : 
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Lemma 2.2 The trajectories associated to 12. 1}) that are contained in the stable manifold and 
satisfy ($, W, \f') P s satisfy that either they are defined for all t € R and \2.2]) holds with 

$(V) 1 

lim -¥ = -■=, (2-6) 
or, alternatively, there exists a t* > — oo smc/i £/iai \2. 3\) holds with 



Proof. That either (|2"2|) or ([23]) hold is just the statement (w,) of Proposition O Then, (|2l)j) 
follows from (|2.2p and integrating the equation. Indeed, for all $ > there exist a r with |to| 
large enough so that for all r < r , then 

4>(t)>$ >0, ^1<-0> <0, ^><f> >0. 

Thus 

d 3 $ 1 

and integrating this expression with t < tq < wc obtain 

1 A {r-pl < Ht) k $(ro) _ ^o) . , (r - ^ (r - r ) 3 



$ 3 y 6 w v uy dr v u ; dr 2 2 6 

for | To | large enough. Then dividing by — r 3 /6 and taking the limit r — > — oo implies ()2.6j) . since 
To can be made arbitrarily negative and <£>o arbitrarily large. 

In order to prove (|2.7p we use the phase-plane analysis of the Appendix 

We employ the transformation (|A.4[) with £ replaced by r for ()2.1j) (see also j7j ) that gives the 
system 

c?$ du 1 o dv 5 , , 

— — =u+-u 2 , — = i + - MU -$ 3 , (2.8 

<22 az i dz i 

which corresponds to (|A.5[) with $ = in the last equation. For further reference, the flow field of 

the phase plane of (|A.5[) is also depicted in Figure Q] where, in particular, the direction of the field, 

the critical point (u e ,v e ) and the attractive separatrix v = v(u) are shown. Using Lemma I A. 3 1 

and the behaviour of trajectories of the system (|A.5[) entering the only critical point (u e ,v e ) as 

z — > — oo it is easy to show that if 

lim $(2:) — > 0, and ||(u,i>)|| is uniformly bounded as z — > —00 , (2-9) 

z— > — 00 

then, by a bootstrap argument, the trajectory (u,v) remains close to (u e ,u e ), in particular the 
estimate (|A.6|) holds for z large enough. That (|2.9j) and that r* = limz^-oc r(z) are satisfied is a 
consequence of the proof of (|2.3j) in [7] and the transformation (|A.4j) . 

Using the first equation in (|2.8p one obtains that there exists positive constants C\ and C2 
such that 

e C2eX *e u ^ z - zo) <S>(z ) < $(z) < e CieXz e u " {z - zo) ^(z ) , (2.10) 

with A = Re(A + ) = (7/2) (1/ (15^ ) > 0, for all z < zq where zq < and \zq\ is sufficiently large. 
And this in particular implies that 

lim $(z)e-"= z = $(z )e-"= Z0 . (2.11) 

Using now the last equation of (|A.4j) and (|2.10p one can infer that, considering r as a function of 

z 

lim (r(z) - T*)ie- U °* = (^) * <Z>(z ) e - u ^> , 

z-t-oo \o4/ 

(where we use u e = (5/9)^ to compute the explicit coefficient). Finally, this and (|2 . 1 1 1) imply 
(EH). ■ 
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2.2 A dynamical systems approach 



As anticipated in the Introduction, in this section we reformulate the main result in terms of a 
Dynamical Systems approach. We first transform (JTTTJ) into a suitable systems of four autonomous 
ODEs, and (jl.2p into its corresponding boundary conditions. 

Since we are interested in solutions for which remains bounded for all £ it is convenient 

to introduce the following change of variables 



H(0 



where the variable r is defined by means of 



1 



« 2 + l)4 



-$(r) 



(e + l)°d£ = dT, r= (r/ 2 

Jo 



dr). 



With this transformation, we have that 
dH 2 £ 

d 2 H 



4 $ + £ 2 + l — 
3 (£ 2 + 1)3 dr 



5 e2 



2 1-K 



rf£ 2 



and (jl.ip becomes 



with 

F(r) 



16 



3 (£2 + 



1 - 



3 (£ 2 + 1)1 
1 CM 2 



d<f> 



9(£ 2 + l)i dr 



..sd 2 ® 

1 I 9 



$3 £( T )2 

14 g 2 
9 £ 2 + l 



--F(r), n€N, se[-l,0] 
1 /208 £ 2 

(£ 2 + i)¥ V"8Tc 2 + i 



10\ d$ 2 
9~ 



(2.12) 

(2.13) 

(2.14) 
(2.15) 

(2.16) 
d 2 $ 



dr 3(£2 + l)¥ dr 2 ' 
(2.17) 

where £ is given as a function of r by means of (|2.13|) . In other words, we use r as independent 
variable, while £ becomes a dependent one, making the system autonomous. 

It is convenient to transform £ further into a new variable that takes values in a compact set, 
namely, we define the variable 9 by 



£ = tan ( 



7T 7T 

2' 2 



Finally, we can reformulate ([2~T3l> and (f2~T6| - (|2~T7| as 



_ 1 

'208 
8T 



— = (cos 0) 



- 1 - (a - l)(cos6>) 2 
(sin*?) 2 



16 224 

-sin^-— (sine) 6 ) (cos0)*$ 
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(cos 0) "5 W + - sin 6>(cos 6) s> * 



that has critical points 



p_ := 1,0,0,-- and p+ := 1,0,0,- 



(2.18) 

(2.19) 
(2.20) 

(2.21) 
(2.22) 

(2.23) 



We aim to prove the following theorem: 
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Theorem 2.3 There exists a heteroclinic connection of the system \2. 19\) - £2.22\) between the 
points p- and p+ given in \2. 23\) . 

We notice that Theorem 1 1.1 1 is just a corollary of Theorem 12. 3) this is implied by (|1.2|) and f)2 . 12[) . 

We point out that the system (|2.19j) - (|2.22j) reduces to (|2.1j) on the subspaces 9 = —tt/2 and 
9 = +7r/2. We shall take advantage of this fact in some of the arguments that follow. 



2.3 Existence of the centre-stable manifold 



We now proceed to describe in detail the construction of a centre-stable manifold at p+ that 
we denote by V+. Let us first define a set of transformations F T (x), r G M. for any given x G 
E+ x K 2 x [-§ , | ] by means of 



(<Z>(T) > W(T),*(r),e(T))=F T (x) 



(2.24) 



where (<S>,W,^,9) solves (l2~T9l - (l2~^2| with ($(0), W(0), *(0), 6(0)) = x. Classical ODE theory 
ascertains that the family of transformations F r (-) is well defined in some suitable interval r G 
(ti(x), t 2 (x)). We have the following result. 

Proposition 2.4 There exists a two-dimensional C 1 manifold V+ contained in the ball Bg(p+) fl 
M 3 x [— tt/2, 7r/2] for some 5 > sufficiently small, tangent to the subspace spanned by the vectors 



(2.25) 





( 3 " § ^ 








-3~* 







Vl = 











I J 




^ 1 / 



// x G V+ , the flow F T defined in \2. 24\ l is defined for any r > and 

F T (x) G V + for any r > 

with 

lim F T (x) = p+ . 



(2.26) 



(2.27) 



Proof. In order to apply standard results it is convenient to extend the range of values of 9, 
replacing cos 8 by | cos 6*| , where the system (|2.19p - (|2.22[) is defined. The resulting system can 
be defined in a neighbourhood of p + and the right hand side of (|2.19l) - (|2.22[) is in x 
(— 7r/2, tt/2)). Since y > 1 we can apply the results in [10]. In our setting, this means the 
existence of a two-dimensional manifold V+ G C^(R 3 x (— 7r/2, tt/2)) tangential to the plane 
spanned by {vi,Vi} at p+ that remains invariant under the flow F T if one can prove that the 
corresponding trajectories on this manifold remain inside a ball B$(p+) for some small 6 > 0. Let 
us show that V+ is invariant . 

Let us consider a four-dimensional cube Q = [1 — 5/2, 1 + 5/2] x [—5/2, 5/2] 2 x [tt/2 — 5/2, tt/2] 
contained in a ball Bs(p+). The cube has four pairs of parallel -3 dimensional- sides. One pair 
with normal direction v±, another pair with normal direction v^, the other two pairs of parallel 
sides contain a plane parallel to the one spanned by v\ and V4. The set Q H V+ gives four C 1 
curves and, due to the tangency of V+ to the plane spanned by v\ and V4, two are contained in 
each of the parallel sides of the cube that are orthogonal to Vi, and the other two are contained 
in parallel subspaces orthogonal to £4. More specifically, one of later is contained in the subspace 
R 3 x {9 = tt/2}. Notice that V+n (R 3 x{9 = tt/2}) gives a portion of the stable manifold associated 
to (j2"TTj) for 5 small enough. Therefore, if x G V+n (M 3 x {9 = tt/2}), F t (x) G B s (p+) for arbitrary 
values of t > 0. On the other hand, for the curve contained in a subspace with constant 8 < tt/2 
and orthogonal to u 4 , we use the fact that 9 is increasing, thus trajectories could only scape the 
cube through the other boundaries that intersect V+. But the points x on the other two boundary 
curves satisfy vi ■ (x —p+) = ic5 for some c > (small or at most of order one). We then use that 



d_ 



(Si • (F r (x) ~ p+)Y' 



(Si • (x- P+ )) { vi 



T = 



dF T (x) 
dr 



T = 
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Since the manifold V+ is tangent to the plane spanned by v% and i>4 it follows, using (Hi) in 
Proposition 12 . 1 1 as well as (|2.25[) that 



I ( ( ^- ( ^ ) ' P+))2 ) | = -3* (n • (x - ,+ )} 2 + o (<* • (* - P + )?) . (2.28) 



Therefore, if 6 is sufhciently small this quantity is negative and the trajectories in V+ remain 
always in the ball Bg(p + ) and (|2.26[) follows. It only remains to show (|2.27|) . To this end, we 
observe that (|2l>2"]) implies lim r ^ 00 6(t) = it/2. Using (j2~28)) we then obtain (|2~27l) . ■ 

For further reference, let us denote by II C M. 4 the affine plane spanned by the stable eigenvec- 
tors at p+, namely, 

n = {w = (u, a) G p+ + K 4 : w - p+ = uvi + av 4 v , a G M} , (2.29) 

with Vi and va as in (|2.25[) . Every w G II can be identified by its coordinates, thus we write 
w = (v, a) and p+ = (0, 0) with this set of coordinates. Since V+ is tangent to II at p + , there exist 
local diffcrentiable parametrisation of V+ . 

Lemma 2.5 (Local parametrisation of V+) Let H be given by V2.29\) . There exists a 8$ > 
and a differentiate mapping A : n -> K 4 that maps a neighbourhood of H into V+nBs o (p+)n{0 < 
tt/2}. Moreover, d v A(0,0), <9 CT A(0,0) G II. 



3 Analysis of the behaviours (11.61) and (11.71) 
3.1 Stability 

We now prove that both asymptotic behaviours (| 1 .6[) and (|1.7[) represent two disjoint open sets of 
solutions of (|2.19|) - (|2.22|) . More precisely, we have the following results: 

Lemma 3.1 Suppose that F T (x) is a solution of \2.19)) - V2.22]) with x G R + x K 2 x (— 7r/2,7r/2). 
Let us also assume that for such a solution lim T „ ) ._ 00 $(r) = oo. Then, there exists a S = 
S(x) > sufficiently small such that for any y G Bg(x) n (R + x R 2 x (— tt/2, 7r/2)) F r (y) = 
($(t),W(t),-&(t),6(t)) satisfies 

lim $(r) = oo . (3.1) 

Proof. It is convenient to use, in order to prove the result, the original equation (jl.l|) that is 
equivalent in the set M + x M 2 x (—tt/2, tt/2) to the system (|2.19|) - (|2.22|) by means of the change 
of variables (pX2|) . (|2T3]) . 

We first recall that ([2~2"2"|) implies that 9 -> -tt/2 as t -> -oo. Therefore, by ([2TT5]) , 
limT-^-oo £ = — oo. On the other hand our hypothesis on $ as well as (|2.12[) and (|2.13[) imply that 

lim |£|iff(0=oo, (3.2) 

^ — y — oo 

then this and (JTTTJ) yields the existence of a £o = £,o{ x ) < with |£ | large enough such that 



d 3 H e 



2 



< -— for all C < Co ■ 



d£ 3 " 2 

Integration this expression gives 



2 - fi fi W£2 



d£ 2 - 6 6 dC 

/ C 4 , Co , ^ 2 #(£o) , ^ , t , dHfa) 
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Therefore, there exist a £1 = £i(x) < with |£i| large enough such that 

We assume, without loss of generality, that £f + a > by taking £1 even larger is necessary. 

Now for F T {y) = ($(r), W(t), tf (t), 0(t)) with y £ B s (x)f] (R+ x R 2 x (-tt/2, tt/2)) we define 
= (£ 2 + l) _1 / 3 l>(r). Since the changes of variables ((2~T2"j) and ((2~T3")) are smooth, it follows, 
using (|3.2[) and (|3.3p and standard continuous dependence arguments for ODEs, that 



Integration of for the unknown /J and p.4[) imply that for all £ < £ 



i 



23 Z/^ 1 Z^ 1 Z"^ 1 Z/^ 1 Z^ 1 Z^ 1 rfsQ 

H(0 > — —j- + / / (4 + a)ds 3 ds 2 d Sl - — 3 — ds 2 dsi . (3.5) 



(£i+a)i ' J s J S1 J S2 Jt J S1 J S2 (H(s 3 )) 

Suppose that 

2* 

H(s)>— -r for (3.6) 

{s z + a) 3 

Therefore, it would follow from (13.51) that: 



23 1 f^ 1 f^ 1 f^ 1 

m)>— 2 7T+o / / (4 + a)ds 3 ds 2 d Sl (3.7) 

(£ 2 +a)3 2^ y si J S2 

where we use that £ 2 > £ 2 . We can then extend the inequality (|3.6I) to a larger range of values 
of £ and therefore the inequality (|3.7j) also follows for all £ < £i with £ in the extended interval. 
Since the integral term on the right-hand side of (|3.7j) tends to infinity as £ — > — oo, we obtain 
(|3~T1) as well. ■ 



Lemma 3.2 Suppose that F T (x) is a solution of (KWjl - (MEjl with x e R+ x R 2 x (-tt/2, tt/2). 
Let its aZso assume that there exists a t* > — oo sucZi Z/iaZ lim r _>( Tji )+ $(t) = 0. Then, there exists 
a 5 = 5{x) > sufficiently small such that for any y G B$(x) H (R + x R 2 x (—tt/2, tt/2)) there 
exists a f * > — oo s«c/i that F T (y) = (^(t),W(t),^(t),6(t)) satisfies 

lim $(t) = . 

Proof. As in the previous proof, it is more convenient to use the original formulation (jl.lj) . We 
again use the smooth transformations (|2.12j) and (|2.13[) to interpret the results between either 
formulation. Thus let H be the solution of (jl.ip associated to F T (x). Let also be defined by 

T * = fo* ( r l 2 + ty^drj, We observe that r* > — oo implies that > — oo, and the hypothesis on $ 
becomes 

lim H(n = . (3.8) 

Thus in regions close to we expect that the solutions are described by and we employ 

the change of variables (cf. Appendix |A1 (|A.1[) and (jA.^ - ). namely, 

^=H~iu, f*=H-lv, H(0 = H(0, £ = fi(z) (3.9) 
where f2(z) is defined by means of 

ttttV- (3 - 10) 

n( z ) (-ff(s)) 3 
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Then, u(z) and v(z) are defined for any z > where z* is given by = Sl(z*). Notice that \z*\ 
may or may not be finite. Moreover, (H, u, v) satisfy 



dH du u 2 dv 5 , , . ,„ _ . 

- r = uH, —=v + — , ^1+ uv-(n 2 + a)H i (3.11) 
Gte dz 6 dz 3 

where all functions, including f2, are functions of z. The hypothesis on $ translates into 

lim ((O(z)) 2 + ft )(i?(z)) 3 = 0. 

*->(*.) + 

The phase-plane analysis associated to (|3.11[) with H(z) = is included in Appendix \M 
Relevant to the current analysis are Lemma [A.2I (where v = v(u) is defined) and Lemma [A. 31 that 
describes the overall flow. 

We claim that there is a sequence {z n } such that z n — > (z*) + as n — > oo and that (u(z n ), v(z n )) £ 
{(u,v) : u > , v < 0} for all n large enough. Before we prove this we note that for any se- 
quence {z n } such that z n — > (z*) + asn-> oo, the trajectory (u(z), v(z)) must be in the half-plane 
{(it, v) : u > 0} for z < z n if n is large enough. Indeed, otherwise the first equation in (|3.11[) 
implies that dH(£)/d£ < for all £ near and this contradicts (|3.8p . Let us now prove that we 
can select such a sequence and that it also satisfies v(z n ) < for all n large enough. 

Let {z n } be such that z n — > (z*) + as n — > oo and suppose that w(z n ) = 0. Then, the third 
equation in p. lip implies that 

^(z n ) = l-((n(z n )f + a)(H(z n )) 3 
dz 

and since the last term converges to zero as n — > oo, it follows that v(z) becomes negative for 
some z < z n close to z n for n large enough. Thus we can construct another sequence {z n } with 
z n < z n , z n — > (z„) + as n — > oo and such that v(z n ) < for n large enough. 

Suppose now that v(z n > for large enough n. Then, the second equation in (|3.11l) implies 
that (u(z),v(z)) arrives to the half-line {u = 0, v > 0} at some z n < z n . For otherwise, the 
last equation in p.lip implies that (u(z), v(z)) crosses the line {v = 0}, and the argument of the 
previous case applies. Therefore, there exists a sequence {z n } with z n — > z* as n — > oo such that 
one of the following possibilities take place: 

lim v(z„) > v{0) , , (3.12) 

n— >oo 

lim v(z n ) < v{0) , , (3.13) 

n— >oo 

lim v(z n ) = v(0) . . (3.14) 

n— > oo 

In the case (|3.12p . we can approximate the evolution of (u(z),v(z)) in intervals of the form 
z G [z„ — L, z n ] by the system (|A.5|) using standard continuous dependence results and Lemma POI 
implies that (u(z),v(z)) enters {(u, v) : u > 0, v < 0} at some z < z n for n large enough, and 
the claim follows. 

Suppose now that (13.131) takes place. Using again continuous dependence we obtain that 
(u(z), v(z)) £ {(it, v) : 1 + 5uv/3 < , u < , i> > 0} = R5 for some z < z n and rt large enough. 
In this region, and with z close to 2*, then u increases for decreasing z. Therefore, d 2 H(£)/d£ 2 
remains positive and dH{^)/d^ is negative as long as (u(z),v(z)) stays in R 5 . Moreover, due to 
the second equation in (|3.1ip \u(z)\ increases for decreasing z. This implies that the inequality 
1 + 5wu/3 < remains valid during all the evolution until z = z*, thus also the inequalities 
d 2 H(£ i )/d£ i 2 > 0, dH(£)/d£ < remain valid. However, this contradicts (|3.8p and (|3 . 1 3[) cannot 
hold. 

It remains to study the case (|3.14p . In this case there exist a small L such that for z £ 
(z n — L, z n ) (u(z), v(z)) remains close to the separatrix v = v(u). On the other hand, (u(z), v(z)) 
must return to {(it, v) : u > 0} infinitely often as z n — > z*. Thus the trajectory must remain close 
to v for z close to z*, or otherwise the trajectory enters R5 giving a contradiction as before, or it 
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enters the region {(u,v) : v < — ^} which contradicts (|3.14p . Then Lemma lA.21 (i) implies that 
d 2 H '(£) / 'dt; 2 > and dH(£)/d£ < remain valid during all the evolution for decreasing z < z n for 
n large enough, and this contradicts (|3.8[) . 

As in the proof of Lemma [31] for F T (y) = (^(t),W(t),^(t),9(t)) with y g B s {x) n (K+ x 
R 2 x (-7r/2,7r/2)) we define H(£) = (£ 2 + l)" 1 / 3 $(r) and the transformed functions \u(z)),v(z)) 
by means of the transformations ()3.9|) and (|3 . 10[) with the obvious changes of notation. 

We then notice that, by continuous dependence of solutions on the initial data, if S > is 
chosen sufficiently small then [u(z),v{z)) enters the region {(u,v) : u > , v < 0} for some z 
close to z* and therefore dH(£)/d£ > and d 2 H(^)/d£ l 2 < for some £ close to with £ > £*. 
We have that (£) is small and d 3 H/d£, 3 > as long as H(£) is small. Integrating this inequality 
for £ < £ we obtain that *^r(£) > 0, *3pr(£) < and £T(£) remains small for £ < £ as long as if is 
defined. Then, vanishes for some £* > -co, so the lemma follows. ■ 

3.2 Characterisation 

We now give necessary conditions for the solutions of to cither satisfy that 

lim H(£) = +oo (3.15) 

or that 

lim H(£) = for some > — oo . (3.16) 

Observe that these behaviours imply (jl.6l) and (|1.7[) respectively, for the corresponding function 
<5(t) given by (|2.12[) and (|2.13|) . We start by giving necessary conditions for (|3.15p . but first we 
need the following auxiliary calculus result. 

Lemma 3.3 Given the polynomials 

•y5 y4 y3 y3 

PAY) = + and P 2 (Y) = , 

n ' 60 12 6 ' 6 ' 

then, they are strictly decreasing and positive for Y < 0. Moreover, if A g M satisfies 1 + 2A > 0, 
then 

Pi(Y) + XP 2 {Y)>^P 1 (Y), for Y<0. (3.17) 

If 1 + 2A < then 

PAY) + 2XP 2 (Y) > --3* max{|l + 2AI*, 1} (3.18) 
5 



for Y < 0, but Pi(F) + 2XP 2 (Y) > ifY < 5/2 - ^25 - 40(1 + 2A)/2. 

Proof. The monotonicity properties of Pi and P 2 are just an elementary calculus exercise. The 

y 5 , r 4 yf. 

60 ' 12 C 6 



inequality (|3.17[) is a consequence of the fact that the polynomial—^- + -yj — c^- = Pi(Y) + (c 



1)P 2 (Y) is non-negative and decreasing if c > 0, in particular 



\p 1 {Y) + XP 2 {Y) = \ 



+ (1 + 2A) — 

60 12 1 ' 6 



is non-negative if (1 + 2A) > 0, thus Q3.17P holds. 

If c < (i.e. 1 + 2A < 0) then P X (Y) + (c - 1)P 2 (Y) < in Y g ((5 - V25- 40c)/2, 0). But 
there the polynomial is larger than or equal than the value of the minimum in Y < 0, namely, 

fl<n + - iwn >(l- - ±) (2 - <4 - *=)»>• > 4(1 - 2c,f 



30 60 15 / v ' ' 15 



and (|3~T8f follows. ■ 

We now give necessary conditions for (|3.15[) to hold. 
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Proposition 3.4 Let us assume that there exists a positive constant c\ — ci(a) > and some 
( eM with 



and 



( y ) 



d > 16(2+ |a|) if |e | 2 >-2a, 



(3.19) 



(3.20) 



such that a solution of IP]) satisfies ((^ ) 2 +l+\a\)(H (£ )) 3 > Ci, dff(£ )M < and d 2 H (£ ) / d(, 2 > 
0. TTien / TQ5)) fto&fe. 

Proof. Integrating (|l.ll) three times for £ < £o we obtain: 

rCo /"Co /-Jo / ^ 

l (Sl+fl) -(^)) 3 , 

Given the polynomials defined in Lemma 13.31 and letting, for every £ < £o 5 



/•Co /"Co /-Cc 

ff(0 > ^(Co) + / / / 



ds^ds2dsi 



(3.21) 



y 



^-lif£o>0 

^ + iif£o<o 



then, we can write 



Co /-Co /-Co 



£ •'81 J S2 



(s 2 + a)ds 3 ds 2 dsi = |£ | £ 



l?or 



(3.22) 



Clearly £/|£o| < -1 if £ < Co and £o < and £/l£o| < 1 if £ < £ and £ > 0, thus in either 
case y < and the polynomials are in the range of values considered in Lemma 13.31 We can now 
distinguish two cases. 

Suppose first that £q < —2a. Then a Gronwall type of argument shows that for any £ < £o 
with £ G [-v^Jof, V^R] then H (0 > 1 > 0. Indeed, as long as > 1 then (|3~2Tj) can be 

estimated from below by 



(3|a| + 1)3 



Pi(Y) 



l£o| 2 



-1 )P2(Y) 



(using (f3T22T) \ Then we can apply (j3~18| with 2A = a/|£ | 2 - 1(< -1), hence 



(3.23) 



Pi(Y) 

using this in (|3.23j) yields 



- 1 P 2 (y) > — 33 max 



H 2 
I6I 5 



(3|a|+l)3 5 



4 3i|a|t 



and (|3.19p implies the claim by a continuity argument. 

Let us assume now that £g > —2a. Using (13.22|) we obtain 



/•Co /-Co /-Co it |5 

J J {st+a)ds 3 ds2ds 1 >^-P 1 (Y) 



by Lemma 13.31 Applying now this inequality to (|3.21[) we obtain the following estimate for 
£o > -2«: 



Co /"Co /-Co 



C >/si Js 2 



(#(* 3 )) : 



-ds 2 ds 1 , £ < £ ■ 



(3.24) 
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Now, we can use a Gronwall type of argument to prove that if c\ satisfies (|3.20[) then (|3.24l) 
implies 

We observe that (|3.25[) holds by hypothesis and that it also holds for £ close to £o by continuity. 
Then, as long as (|3.25[) is satisfied, f|3.24[) implies that 



g (fl>( , « , J t Mp, |r )- 'itfftN) affi , f) . (3 . 26) 



|&| 2 + i + M7 2 " ' ci 

We can apply Lemma T3.31 and this implies that the last term in (|3.26|) can be estimated by 
the previous one for any £ < £o if c\ , a and £o satisfy 

d >8(|Co| 2 + l + |a|)/|^o| 2 (3.27) 

and (|3.25p follows. Let us then prove (|3.27j) . 

If | £o| > 1, ([3~2U|) implies ([3127]) , If |£ | < 1 we consider two further cases. For |£ - £o| < 2 we 
obtain that the last two terms in (|3.26p can be bounded from below by —8(2 + |a|)/(3ci). But 
this quantity can be absorbed by the first term if c\ > (40) * (2 + |a|), which is satisfied if (|3.20p is 
satisfied. Therefore p.25[) holds for this range of values. 

On the other hand, if |£ — £o| > 2, then the second term of (|3.26[) can be estimated from below 
by (24/15)|£ - £o| 3 , while the last term in (|3~26| can be estimated by -(2 + |a|)|£ - £ | 3 /(3ci). 
Then, we can absorb the last term in (|3.26p into the second one if c\ > 5(2 + |a| ) / 24 which is 
guaranteed by (|3.20p . 

Thus, the inequality (|3.26p holds for arbitrary values of £ < £o> an d this implies f|3 . 1 5[) by 
taking the limit £ — > — oo. ■ 

We end this section by giving necessary conditions for (|3.16p to hold 

Proposition 3.5 Let us assume that there exist positive constants c 2 and C3, depending on a, and 
some (o£K with 

1 

^-(|£ | 2 + 1 + H)I <1 (3.28) 
C3 iu 

such that a solution of (XT]) satisfies < ((£o) 2 + l + M)(#(£o)) 3 < c 2 , (\Co\ + l + \a\)idH(^ )/d4 > 
C3 and d 2 H{^)/d^ 2 < 0. Then there exists £* € (— 00, £0) such that \3.16}) holds. 

Proof. Suppose that c 2 is sufficiently small. Then, as long as < ((£) 2 + 1 + |a|)(i?(£)) 3 < 2c2 we 
obtain from ([l.ip that d 3 H (£) / d£ 3 > 0. Integrating this equation over (£,£0) once we obtain that, 
as long as (£ 2 + 1 + \a\)(H (£)) 3 < 2c 2 is satisfied for £ < £ , then d 2 H(£)/d£ 2 < and, integrating 
a second time, also that dH(£)/dt; > C3 ( | ^0 1 2 + 1 + |a[) — 3 . Then this concavity implies that i?(£) 
vanishes at some £ = £*. But a third integration implies that 

H(0 < - ^> (& - < fa* - 1 C3($ °' e) 

1 

thus £* >£o~"j|-(|£o| 2 + l + M) 3 - Finally the condition (|3.28p implies that we can replace £0 by 
£ £ (£*,£o), thus (£ 2 + 1 + |a|)(i?(£)) 3 < 2c 2 follows in this interval and the result follows by a 
classical continuation argument. ■ 




4 Shooting argument 

In this section we apply a standard shooting argument to prove the existence of solutions of 
([2~T9|) - ([2~22|) such that (TO) holds, and such that $ remains positive and bounded for all r € K. 
Specifically, the main result of this section is: 
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Proposition 4.1 There exists a solution of \2.19\l - \2.22\) ($(r), W(t), *(t), 0(t)) defined for all 
t £ (-00,00) such that lim r _ >00 ((f>(T), W(t), ^(t), 0(t)) = (1,0,0, 7r/2) and satisfying $(r) > 
for all t £ M, JO)) and 

lim inf $(t) < 00 . (4.1) 

t — y — 00 

The proof of Proposition 14.11 is divided in several steps. First we prove that points placed in 
the curve V+ H M 3 x {6 = -| — e} with e > sufficiently small, yield solutions of the equation p. II) 
satisfying the hypotheses of Proposition [3]4] if v > and those of Proposition 13.51 if v < 0: 



Lemma 4.2 Let So and A(i/,a) be as in Lemma \2.5[ Then there exist v$ > and e > 0, such 
that for w — (y, a) £ II itraf/i i>q < v < So/ A and a = —s the trajectory associated to \2.19)) - H2.22\ ) 
starting at A(i/, a) satisfies il.6\) . Moreover, if —So/4 <v< —Vq and a = —e the corresponding 
trajectory of (KB\) - (2~l2]) satisfies (Tl\) . 

Proof. The dynamics induced by the system (|2.19j) - (|2.22j) on the invariant subspacc R 3 x {9 = 
7r/2} have been summarised in Proposition 12.11 and Lemma 12.21 In particular, the trajectory 
starting at A(i/, 0) with v > sufficiently small satisfies (|2.6[) and, as it can be easily deduced, also 
that 

lim ^ = -i and lim 

r— > — 00 T 2 t— >— 00 T 

Then, classical continuous dependence results for ODEs imply that for any po > arbitrarily 
small and vo > small enough there exists e sufficiently small such that, for vo < v < So /A the 
trajectory starting at A(v, — e) at r = satisfies: 



< Po|to| 3 ■ 



W{r ) + T l 



<Poko| 2 , \^(t )+t \ < Po \t \ (4.2) 



-H(Co) > ci|4o| 3 , — ^ — <0, — ^ — >0 



for some r < 0. Using (|4~2"j) and (|2~T2]) - (|2~T5l) to get £T, diJ/d£ and d 2 H/d^ 2 at the value £ 
(given by (|2.13[1 ). we obtain 

d H(to) d 2 H(£ ) 

d4 ' d£ 2 

where C! > can be made arbitrarily large choosing e sufficiently small and |to| sufficiently large 
to guarantee that p. 201) is satisfied. Then we apply Proposition 13.41 to obtain (|3.15[) and hence 
(|1~d) follows. 

On the other hand the trajectories starting at K(v, 0) with v < satisfy lim r _ J>r + $(r) = 0, for 
some r+ > —00. Moreover, (|2.7p is satisfied, as well as 

1 1 
lu^(r-r^W(r) = l(^y and ^lim (r - r.)«*(r) = ~ 

Suppose now that — <5o/4 < C < — ^Oj 17 = — £• Assuming again that e is sufficiently small we 
obtain that the numbers 

$(r ) /64\3 1 3/64\* « 3/64^ 

(^T)J-Uj ' (-o-n)^(r )- 3 (-j , ( T0 -r^(r ) + -(- 

can be made arbitrarily small for tq close to t* , to > r* . We can use this approximation to obtain 
that 

3 ' de 

where c 2 oc (r — r*) 1 / 4 and c 3 oc (t — ) -1 / 4 , thus they can be chosen to satisfy (|3.28p by taking 
e > sufficiently small. We can now apply Proposition ^. 5l to conclude the proof of the result. ■ 
Next we prove that if for every compact set K C (—00, 00) we have that Uminf T _ v ( Ti \+ $(t) = 
for some r» > —00, then lim r _>( Tjii )+ $(t) = 0. Therefore, we will be in the situation stated in 
Lemma 13.21 and it will be possible to prove continuity of this behaviour for small changes of the 
initial values. 



Hfa)<dM-i, ^>c 3 , ^^<0 
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Lemma 4.3 Let 3>(t) be a solution of \2.16)) - {2.1T^ defined in some interval (r*,r*) with r* > 
-co, t* < oo, &(t) > for t > TTien 

lim inf $(r) = . (4.3) 

r— >(r») + 

implies \1. 7| ) /or i/iis vaZwe r = r* . 

Proof. It is easier to work with the original equation observe that then, (|4.3[) is equivalent 

to 

lim inf = (4.4) 

for 6 given by r* = J Q * {rj 2 + l)^dr/. Let us then prove that (|4.4[) implies ()3.16j) . and therefore 
(fl~7l) will follow. 

We argue by contradiction. We then assume that (|4.4j) is satisfied, but (|3.16|l docs not hold, 
this means that also 

lim sup H(£) > . (4.5) 

On the one hand (|4.4[) gives the existence a decreasing sequence {6} such that 6 — >• as n — >■ oo, 
H(£ n +i) < H(£ n ) and limn^oo i?(6i) = 0. And (|4.5[) implies the existence of a sequence with 
elements 6 € such that H(£ n ) = H(£ n ). Then there exists another sequence {6} with 

6i £ (6n6i) and lim n _ ) . 00 6 = 6j where local minima are attained, i.e. satisfying 

H(U= min #(£), lim (£ n ) = 0, ^^ = 0, > . (4.6) 

Let also {6} be the sequence where local maxima are attained, such that 6+1 < 6 < 6 and 
satisfying 



#(U = max _ ,#(£), lim sup ff(f n )>0, — (6) = 0, — (£„)<0. (4.7) 



dff,« _ d 2 ^ 

£e(£„ +I ,f„) n^cx) at 

Let us now show that 

d 2 #(6) 



lim 



d^ 2 



> . (4.8) 



Indeed, from (|1.1[) we obtain that d 3 H/d^ 3 > —C\ if £ € [6+1,6] and integrating this inequality, 
we also obtain 

ET/V "\ \ uY<* \ 1 ^ H(^ n ) (61 — 60 (61 — 60 

ff(£„) > ff(£„) + — ^ ^ . 

Now, if (gU) fails, it follows that #(60 > for some subsequence, and this contradicts (|4.6I) 

and g7l). Thus (OJ holds. 



We now claim that (|4.8p implies that -ff(£) vanishes for some £ £ [61+1,6]- Indeed, since 
d 3 H/d£ 3 > —C\ we then have that, for n large enough, 

d 2 ff(Q d 2 g(6Q f : l gg^Q , 

This implies that dH(£)/d£ > for n large enough with £ G [6+1:6]: but this contradicts the 
definition of 61+1 > an d so for n large enough there is a first value £ £ [6+ij 6] such that H(£) = 0, 
i.e. H535| holds. ■ 

Remark 4.4 Notice that a classical Gronwall argument implies that any solution 3>(t) can oe 
extended for arbitrary negative values of t as long as $(t) remains away from zero. More precisely, 
i/liminf T ^ T + $(t) > for any tq > r* > —00, it is possible to extend $(t) as a solution of V2.19\) - 
112.22]) for times r > t* — (5 and some <5 > 0. Reciprocally, the maximal existence time, due to 
Lemma \4~^l * s finite and it is given by r* > —00 if lim inf T _j. T) , $(r) = 0. 
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We are now ready to prove Proposition ^. II 
Proof of Proposition |47H We consider the one-dimensional family of solutions of (|2.19[) - (|2.22p 
obtained choosing in Lemma 14.21 the parameters a = e > with e > small enough and v G 
(— <5o/4, <5o/4). We define as U + the set of values of v such that the corresponding solution of 
(I2.19p ~ (|2.22[) satisfies (|1.6|) . On the other hand, we denote by U- the set of values of v such that 
the corresponding solution of (|2.19j) - (|2.22|) satisfy (|1.7|) for some r* > — oo. Due to Lemma |4~2"1 we 
have that U+ ^ and W_ ^ 0. Moreover, by definition U+ P\U- = 0. Due to lemmas l3~Tl and [3721 
we have that the sets U+ and U- are open sets. Therefore, there exists v G ( — such that 
v i U+UU-. 

The corresponding solution of f|2. 19[) - (|2.22[) associated to the parameter v has the property 
that, for any tq > — oo we have inf re ( ro 00 ) $(r) > C_(ro) > 0, since otherwise v G U- due to 
Lemma l4~3l This implies also that sup Te / T0jO0 ) $(r) < C+(tq) < oo because the right-hand side 
of (|2.19|) - ([2.22|) is bounded in compact sets if $(r) > C-(to). Therefore, this solution is globally 
defined for r G (—00,00). Moreover, (|4. 1[) holds, since otherwise v G U+ and the result follows. ■ 



5 Oscillatory solutions 



We recall that the final aim is to prove that the solutions found in Proposition 14.11 have no 
alternative but to approach the invariant subspacc 9 = —n/2 as r — > —00 and they remain 
uniformly bounded while $ stays positive (sec (|1.8I) - (|1.9|0 . The argument is by contradiction and 
in this section we prove the following lemma that is the first step in the argument. 



Proposition 5.1 Suppose that ($(t),W(t),^(t),9(t)) is a solution of H2.19\) - l[2^) defined for 
all t G (—00, 00) and satisfying 

lim inf $(r) < 00 (5.1) 



and that 



lim sup &(t) 



d$(r) 



dr 



d 2 $(r) 



dr 2 



Then, there exists a decreasing sequence {r* } with lim Il _ i . 00 t* 
e„ > small enough such that 



HO 



max 

t£[t*- e n ,T*-\-e T 



Mr) 



id that 



lim sup $(t*) = 00 . 



= 00. (5.2) 

-00 and a sequence {e n } with 

(5.3) 
(5.4) 



Before we prove this result prove three auxiliary lemmas. First we show that there exists a 
decreasing sequence of local minima attained at certain r = t„ with lim Il _ i . 00 r„ = —00. 

Lemma 5.2 Let ($(t), W(t), ^(t), 9{t)) satisfy the assumptions of Provosition [57l\ Then, there 
exists a decreasing sequence {t„} such that lim„_ i . 00 t„ = —00 and that 



lim $(t„) < 1 
Proof. First, we claim that 



and 



dr 



0, 



dT 2 



> for all n . 



lim inf <f>(r) < 1 , lim sup <f>(r) > 1 



(5.5) 



(5.6) 



Indeed, suppose first that lim inf r _ i ._ 00 $(r) > 1. Then, there exists £0 > and To sufficiently 
negative, such that $(r) > l + 2e for r < r . Then, (j2~T2]) and (|2~T3)) imply H (£) > (1 + £ ) |^| ~ 2/3 
for ^ < ^0, where £0 is related to tq by means of (|2.13|) . This inequality applied to (jl.lj) gives 



d 3 H 
7i^" 



^3-(e 2 + a)<-e^ 



(5.7) 
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for £ < £o and some E\ > (by taking a more negative To if necessary). Integrating (15.71) 
three times for £ < £o gives > £i|£o| 5 -Pi(l + £/l£o|) (where Pi is as in Lemma [3.3p . Thus 

lim^_ ! ._ 00 H(£) = oo, but this contradicts (|5.ip . 

We now prove the second inequality in (|5.6[) . Suppose on the contrary that lim sup r _ > _ 00 $(t) < 
1, then $(r) < 1 — £o for some £o > and for r < to if To < with |ro| large enough. The 
transformation (|2.13|) - ()2.12|) . with the obvious correspondence in notation, implies that i?(£) < 
(1 — £o)|£|~ 5 for £ < £qi whence (jl.ip yields 

^ = ia-^ 2 +«)>^ 2 (5-8) 

for some £ 2 > and £ < £ - Integrating (JSTHJl for £ < £o we obtain < -£i|£ | 5 -Pi(l + £/|£o|)- 

This implies the existence of a > — oo such that lim^_ j .^*- ) + H(£) = 0. This contradicts the 
assumption that £* = — oo (see Remark |4~4| and ()5.6[) follows. 
Suppose now that 

lim inf $(r) = lim sup $(t) = lim $(r) = 1 . (5.9) 

r— » — oo r— V — oo r — *~ 00 

We define a sequence of functions {<I> rl (s)} with s £ [—1, 0] as follows. For every n € N the variable 
£„(s) is given by, cf. (|2T2]l . 

(1 + r/ 2 )*d77 , neN, (5.10) 
then, each $ n (s) is defined by, cf. (|2.13|) . 

<f>n(s) = {l + \Us)\ 2 V H(Us)), s<0, neN. 

We observe that then 

$n(s) = $(s-s n ), where s„ = / (1 + r^cfr? , 

J— n 

where $ solves (|2.16p - (|2.17[) . Also, for every s G [—1,0] the corresponding sequence t„ = s — s n 
converges to — oo as n — > oo, since lim n _ i . 00 (s n ) = oo. On the other hand, the functions $n(s) 
solve (cf. (l2T6l) - (l2T7in 

^-^-frf -'•(■)■ "i-i,"! (^-") 

where F n (s) is given by the expression of F in (|2.17p with $ and £ replaced by and £„, 
respectively. 

Consider now the result of integrating (|5.1ip : 

°'°'° r 1 is-wr+'i***, (5,2) 



Si 

r r0 



($„( S3 )) 3 |^( S3 )| 2 + 1 

F n (s 3 )ds 3 ds 2 dsi . 



Si •! S 2 



We now pass to the limit in the integral terms. Observe that the assumption (|5.9p implies that 
linin-^oo 3>n(s) = 1 uniformly on [—1,0]. Moreover, (|5.10p yields lim„_ > . 00 £, n (s) = — oo uniformly 
on [—1,0]. The first term in (|5.12p can be seen to converge to zero using the limit properties of 
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<J> n (s) and £n(s). In the last term we integrate by parts where necessary in order to get integrands 
with $„(s) as a coefficient (this gives boundary terms with a double or single integral, but these 
are estimated similarly, because s G [—1,0]). The resulting integrands have <f> n (s) multiplied 
by a function of £ n (s) and its derivatives, which can be computed using (|5.10[) : d^ n (s)/ds 



(|£„(2)| 2 + I)" 479 and d 2 Us)/ds 



-4(|£„(2)| 2 + l)- 17 / 9 /9. Then, one can conclude that the 



limit of the last term in (|5.12j) tends also to zero as n — > oo, and we are left with 



lim 

n— »oc 



d$„(0) 



ds 



d 2 $ n (0) 



ds 2 



0. 



but this contradicts (|5.2l and (|5.9p cannot hold. Then liminf T _ ) ._ 00 $(r) < limsup r _ i ._ 00 $(r). 
We can now construct a sequence that satisfies (|5.5p . We first take the following quantity 



lim sup ^(t) + lim inf $(t) 



T— ¥ — OO 



(that might be infinite if limsup T _ > ._ 00 $(r) = oo). Due to the continuity of $, there exist 
decreasing sequences {•?„} and {t„} such that limn^oo f„ = lim n ^oo f„ = — oo, f„ < f„ and 
that maxTgjf^^) &{ T ) < We define another sequence {t„} by 

$(t„) = min $(r) , 

and this one satisfies (|5.5|) . ■ 

We continue with another consequence of assuming (|5.2|) . namely, 

Lemma 5.3 Suppose that (^(t), W(t), 'J'(t), #(t)) satisfies the assumptions of Proposition [57l\ 
Then, at least one of the following identities holds: 

lim inf $(t) = 0, lim sup $(t) = oo . 

T— > — OO T — ^ — OO 

Proof. We argue by contradiction. Suppose that 

lim inf $(r) > and lim sup $(t) < oo , 

t— > — oo t— > — oo 

then, there exists a Co > and a to sufficiently negative such that 

-^<$(r)<C , t<t . 
Oo 

Then, (|5.2[) implies the existence of a sequence {f„} with f„ — > — oo and such that 



(5.13) 



lim 



d$(f„) 



r/r 



d 2 <J>(r n ) 



dr 2 



(5.14) 



We now define a sequence of functions $„(2;) by means of <& n (z) = §{z + f„) and observe that 
they solve (|2.16[) - (|2.17p with the obvious changes in notation and with £„(z) defined by 



z + r„ = 

They also satisfy, due to (|5.14j) . that 

lim 

This allows us to introduce, for every n, the length scale 



(1 + rj 2 )^ drj , 77, eN. 



d$„(0) 



rf 2 $n(0) 



dz 2 



7>i = 



<P* n (0) 

dz 2 
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which clearly satisfies lim„_i. 0O -f n = 0. We now set z = jnZ, £ n (z) = £, n (jnz), $n{ z ) 
and F n (z) = F n (j n z), to obtain that i>„ satisfies 



d 3 $ n 

dz 3 



(7n 



1 



s 2 



£ 2 + 1 



+ { ln fF n {z), neN. 



It is clear that there exists a C > such that for n large enough 

d$ n (0) d 2 $„(0)\ = fd* n (0) d 2 <§„(0) 

1 < O and 



*«(0),- 



f/2 



dz 2 



\ dz 



dz 2 



> 



C 



(5.15) 



We can now use classical continuous dependence results for ODEs. Let $oo denote a function that 
solves the limiting problem d 3 $ oc /dz 3 = with initial conditions close to ($ n (0), d$ n (0)/dz, d 2 $„(0)/d 
for n large enough, and thus also satisfying (|5 . 1 5[) (with n = 00), then 

$„(z) -> $oo(z) as ti ^ 00 , 

as well as its derivatives, uniformly in compact sets of z. 

Since $oo is a polynomial at most of second order that is not identically constant, there exist 
values of s such that either $oo (z) = or & x (z) > 2M. But this contradicts (|5.13|) . ■ 

The third auxiliary lemma is the following: 



Lemma 5.4 Suppose that ($(t), W(t), 'I'(t), #(t)) satisfies is the assumptions o f Provosition [377 
Then, 

lim sup $(r) = 00 . 



Proof. Suppose that 



lim sup $(t) < 00 . 



(5.16) 



T— f — OO 



Then, due to Lemma 15.31 we have liminf T _ ) ._ 00 $(t) = 0, and there exists a sequence of points {t„} 
such that (|5.5p and that limn^oo $(t„) = 0. For every n we introduce the changes of variables 



$(r n ) , s = e„ 3 (r - r„) and $(t) = e n ip(s) . 



(5.17) 



Observe that if £ n is related to r„ by means of (|2.13|) . In particular, this reflects that when |r — t„ 
remains bounded as n — > cxj, then |£| — > 00 as n — > 00, a fact that we shall apply below. 
We write (|2.16|) in the new variables, then <p solves 



with (cf. 12171) 

F n (s) 

and subject to 



16 



dV 

ds 3 



i 2 



a 



T 



3 (£ 2 + l)* 



'"^ + 1 

t2 



1 



14 £ 2 



9 



(£ 2 + 1)* 



208 i 2 



p(0) = l, 



81 £ 2 + 1 

d(^(0) _ 
ds 



f 1 
10 
~9 



&n F'n (#) 

<p(s) 

dip(s) '. 



(5.18) 



d 2 ^(s) 



0. 



ds 



d 2 y(0) 
ds 2 



3(^2 + 1 )f d s 2 



> 0. 



(5.19) 



The coefficients involving £ are functions of s; £(r) = £(t„ + Sils). Then, since e n — > 0, we can use 
classical continuous dependence results for ODEs to approximate the solutions of (|5.18[) - ()5.19l) by 
the solutions (p of the limiting problem (with e n = 0) 



d^_^ = (} 

ds 3 tp 3 



<p(0) 



d<p(0) 
ds 



d 2 <p(0) <jP<p(0) 



ds 2 



ds 2 



> 
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on compact sets of s. Clearly, d 2 (p(s)/ds 2 is increasing for s > 0, and there exists a S > such that 
d 2 (p(s) / ds 2 > S > for all s > 1. Thus, d<p(s)/ds increases at least linearly for s large enough and 
one can find a value sq > such that d(p{sa)/ds > 2. Continuous dependence results then imply 
that d<p(so)/ds > 1 if n is sufficiently large (see (|5.18[0 . We now return to the original variables 
and get estimates on $(r) for every n in an interval around the local minimum. Using (|5.17p . 
(|5.16|) and the notation 



we obtain that 



lim <&(t„) = , lim 



dr 



•En s : 



oo , 



,. d 2 $(f„) 

lira — — = co . 

n-¥oo dr 



For n large enough and r such that |r — t„\ < 1, there exists a positive sequence {B„} with 
lim^-^oo B n = oo and such that 



B n < 



1 



$(t)3 £(r) 2 + l 



< B. 



n+l 



Let us denote by 



G„(r) = ^ jf 1 jf 2 F(T 3 )dT 3 dT 2 d Tl , 
where the function F(r) is given by (|2.17p . and let also 

Pn (r) = *(r B ) + «^(r- r n ) + 
Then, integrating f|2 . we can write for re [r n — 1, T n + 1] that 



B n d k (r - f„) 3 < d k ($- Pn -G n )(T) < B n+1 <* 



fc = 0,1,2. 



3 dr fc dr fe 3 dr" 

Then, taking n large enough and combining the inequalities (|5.20[) we obtain that 



(5.20) 



|d>(r)| <d>(f„) + 2^^|r 



dr 



d$(r) 



dr 



dr 2 

f'2 f 



„d$(f„) „d 2 $(f n ), 
< 2 — + 2 — -V^ r - f n 



+ |G„(r)|, 



dG„ 



dr 



d 2 $(r) 



dr 2 



dr 2 



< 2 



d 2 $(f») 
dr 2 



dr 
d 2 G n 



(5.21) 



dr 2 



and that 



$(r)>$(r„) + ^^(r-r„ 



d 2 $(f») (r - f ra ) 2 
dr 2 2 



+ B rl (r-f„) 3 + G„(r) 



(5.22) 



for some sequence B n > with lirrin^oo B n = oo and for all r G [r n — 1, r„ + 1]. We now observe 
that F has the form F(r) = /i(r)$(r) + f 2 (T)d&{T) / dr + f 3 {T)d 2 ^{T) / dr 2 where the functions 
/i( r )i /2( T ) and /a(r) converge uniformly to zero on sets |r — f„| < 1 for every n (cf. (|2.13p and 
(|2.17p ). This allows to get estimates on the integral terms (those involving G„) as follows: 



\G n (r)\ 



dG„(r) 




d 2 G„(r) 




dr 




dr 2 





--f„l<l 



d$(r) 



dr 



sup 

|t-t„|<1 



d 2 $(r) 



dr 2 



for all n with r such that |r — r„| < 1 and where i n — > + . Applying this to (|5.2ip and to (|5.22p 
we obtain, taking n sufficiently large, that 

*(r) > d ^A^± + *^(t-t b )» + ^ _ ^ 



dr 2 
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for |r — f„| < 1. Choosing, say r — r„ = 1, we obtain that liuin^rx, $(f n + 1) = oo, but this 
contradicts (|5.16|) . whence the lemma follows. ■ 

We are now ready to prove Proposition 15. II 
Proof of Proposition ^, ll The assumptions imply that we can use the statements of lemmas l5.2l 
and 15.41 In particular, by Rolle's theorem, we can guarantee the existence of local maxima in each 
interval (r n +i, T n ) where {r n } is the sequence of minima defined in Lemma 15.21 We then observe 
that the regularity of a solution $ of (|2.16j) guarantees that the points at which $ attains local 
maxima or minima are isolated. Otherwise $ would take constant values on closed intervals, but 
constants are not solutions of (|2.16[) . Hence, we can define the sequence such that (|5.3p holds. On 
the other hand, Lemma T5.4I implies (|5 .4[) . ■ 



6 Properties of oscillatory solutions 

In this section we proof the following proposition: 

Proposition 6.1 Let the assumptions of Provosition [57l\ hold and let {t*} be the sequence found 
in this proposition. Then there exists uq G N large such that for all n > no 

(6.1) 

We observe that this result is in contradiction with (|5.4[) . We now define a sequence {r™ m } as 
follows 

Jfot min\ J2 ((, I rain \ 

Hrr n )= min $(r) , V ' =0, ) \ > >0, (6.2) 



i.e. $ reaches the minimum in the interval (r*, t*_^ at t = r™ 1 '" 1 . Then as part of the construction 
necessary to prove (|6.1[) it will follow that: 



Proposition 6.2 Let the assumptions of Provosition HOI hold. Then the sequence {t™™} given 
in is well-defined, limn^oo r T """ = — oo and there exists no such that for all n > no then 

d>(C!?) < $(r™ n ) . (6.3) 

The proofs of these propositions are divided in several steps that we outline below for clarity. 
We first identify a two-parameter family of polynomials that approximate H near a large maximum 
of $. Most part of this analysis is done in Appendix [B] where we identify and give some properties 
of the polynomials that solve (|1.10[) (with the reverse sign). Then for each n and around r* we 
identify a length scale that transform these polynomials into polynomials of order one. We then 
translate the properties found into the rescaled polynomials. We also rescale accordingly the 
function H near each £* defined by (|2.12j) for r = r* and give the approximating lemma that 
in particular implies that H will get close to in a linear decreasing way. Next we adapt the 
matching lemma, Lemma IC.ll in Appendix \C\ that gives the behaviour of the solutions in the 
inner regions near each £™ m = £( T ™ m )- From this result we can conclude that the approximating 
polynomial in the outer region must have a double zero in order to match. This, in particular, 
reduces the class of approximating polynomials to a one-parameter family. We finally derive an 
iterative relation between the elements of the sequence {$(r*)} if n is large enough that implies 
Proposition 16. 11 as well as a relation for the elements of {^(t^"")} that implies Proposition ^. 21 

6.1 The outer variables and the auxiliary polynomials 

Given the sequence of {t*} found in Proposition [5T] see (|5.3[) . and the sequence of local minima 
{t™ 11 } defined in (|6.2[) . we define the sequences {^}, {CT m }i {^n} an d by means of (see 
(1531) and jEl): 

(rf + l^dn, r™ n = / {n 2 + l)*dn, (6.4) 
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iapf — ( } 

Observe that the definition of r* implies that f3 n < 0, that lim n _ i>00 £* = — oo and that linin^oo £,™ m ■ 
— oo. 

Observe that M n is the value of the maximum of $ at each r* rescaled appropriately with 
the position of the maximum in the variable £, £* (this scaling near a maximum resembles that 
<f> ~ there and that H is approximated by a fifth order polynomial, cf. (|B.3|) and (|B.4I) 

in Appendix [B]) . The definition of j3 n results from similar considerations, but is a parameter that 
captures the value of the second derivative of <J> at r* . 

Following Appendix [B] associated to every t* we construct the two-parameter family of poly- 
nomials 

P(Z n ;M n ,f3 n ) with Z n =± 

that solve (|B.3[) with Z replaced by Z n , and that are given by (|B.6|) . We recall that they satisfy 
that P(l; M n , p n ) = M n . We shall see later that these polynomials are close to — (Ci) 5 -^(£) f° r £ 
close to £* in an interval contained in [£™i,£™ m ]- Thus we have to consider Z„ in some interval 
containing Z n = 1 and where P stays positive. Moreover, since £* — > — oo the approximation will 
be applicable for Z n > only. In this regard, for each n, we have derived a number of properties 
that arc outlined in lemmas TB. 11 IB.2[|B.3l and lB.4l These give, in particular, that the largest root 
of P(Z n ; M n , f3 n ) in Z n < 1 is attained at a value Z n = Z${M n , /3„) for every M n > and every 
fin < 0. It is also shown that for every M„ > there exists a unique value /3 n = /3*(M„) such that 
P(Z n ; M n , /3*(M n )) has a double zero at some if„ = Z*(M n ) > 1. 

In these lemmas the asymptotic behaviour as M„ — > and as M n — > oo of j3*(M n ), Z*(M n ) 
and Zo(M n ) = Zo(M n , /3*(M n )) is also given. But, as we shall see later and assume now, the 
sequence {M n } is bounded. 

Taking these considerations into account, we now introduce a rescaling of P(Z n ; M„, f3 n ) for 
every n in order to have values of order one in the relevant range of parameters. Namely, we set 

P(U; M n , AO = r^ fM , Cn = _J_ ( 1 _ i) . (e.e) 

1 /3 

Observe that now the variable £ n is meaningful in an interval around £ n = and with £„ < 1 /M n ■ 

i 

We note that M„ is a characteristic length scale which measures the distance between Z n = 1 and 
Z„ = Z„(M„), relevant if M n is very small, see Lemma TB. 2 1 We also observe that the polynomials 
P(C„; M n , fi n ) are explicitly given by 



P(C«; Mn, /3„) = - jg (mIc' - 5M*Cn + 10) + Qm| C + |m|c« + l) + A: y ■ (6-7) 
It is natural to define the following values of £ n : 

Co(M n ,fi n ) = - Zo(Mn ^ 1 n) - 1 , (6.8) 

(see Lemma EH (|B~14} ), thus clearly ( (M n ,(3 n ) = min{C„ > : P(( n ;M n ,p n ) = 0}. And if 
fin = fi*(M n ) we define also 

UMn ) = - Z ^ Mn }-\ (6.9) 

therefore P has a double zero at this value (see Lemma fB. II) . When fi n = fi*(M n ) and for simplicity 
of notation, we shall write: 

P(C„; Mn) = P(Cn, M n , fi*{M n )) 

and 

Co(Af n ) = Co(M„,/3*(M„)). 
The following result follows easily: 
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Lemma 6.3 For every n and M n the polynomial P(Cn', M n ,/3 n ) solves 

^%^ + (1-m| C ) 2 =0 (6.10) 
etc, 

with initial conditions 

p, n ^ \ -j dP(0;M n ,*3 n ) _ mk d 2 P(0;M n ,p n ) _ p n IOmJ 

We now reformulate the results of Appendix iBl for these approximating functions: 



Lemma 6.4 Let P(Cnl M n , (3 n ) be given by |fi.7| ). They satisfy that if ' j3 n > (5*(M n ), then -P(C™; M n , j3 n ) > 
in d < and if /3 n < j3*(M n ) then there are two zeros of P(C, n ; M n /3 n ) in £„ < 0. The derivative 
°f P{Cn', M n , j3 n ) with respect to £„ is positive at the largest root in £„ < 0. If f3 n — /3„(M„) there 
is only one double zero in £ n < and is placed at £„ = ( t (M„). Moreover, 

C*(M n )~ -12^(1 + m|), &(M n )~-(^V as M n ^0 (6.12) 



V 2 

d 2 P(UM n );M n ) f3\ 

as M n — > . (6.13) 



ac 2 V2 

Also, the value i6.8\) is well-defined for every n and M n , and if fj n = /3*(M n ), then 

«M„>~(i) J , «i«.-(|) S „ M„.0. (,», 

Finally, we have also that 

Lemma 6.5 T7ie waZwe Co(M n , /3 n ) is the unique root of P(Cn] M n , /3 n ) in £ n > 0. Moreover, if 
i 

Co(M n , f3„) < 2/Mn , there exists a positive constant cq independent of M n and (3 n such that 

dP(Co(M n , /?„); M n , f3 n ) 



d( 



< -c max{l,M r ?} . 



Remark 6.6 We point out that the case Mq — > corresponds to the the approximating polynomials 
obtained for jl.5\) in J^j. The asymptotics \6.1J$ are in agreement with this observation. 

6.2 The sequence of rescaled H(£) near each £* 

In order to compare H with a polynomial P(Cn', M n , f3 n ) we need to apply the scaling (16. 6[) to H 
around £ = £*. We then obtain: 

Lemma 6.7 Let us assume that $ satisfies the assumptions of Proposition \5A\ so that the se- 
quence i f 5. 3\) is well-defined. Let H(£) be the solution of 11. 1]) related to $ by means of &2.12]) and 
\2.1S\) . Let the sequence of functions {'Hn(Cn) i Cn S R} be defined by 

** |5 jt ir n t ft- \ ^ I ^ 



H(0 = |£| 5 M„W n (C„) , Cn = r £ - 1 . ( 6 - 15 ) 

Ml / 

i/ien, /or eac/i n, H n solves 

^ +ft , (( „ ;M „) = ^, (6.16) 

w/iere 



, i?n(Cn; M„) = (1 - MlCnY + . ( 6 - 17 ) 
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with initial conditions 



n n (o)= 'S'', ^M = |m| ^ 4 (6.i8) 

(l+lel 2 )* rf Cn 3 (|a| 2 + l)^ 



d 



rf 2 H„(0) _ /3n |g|j ^ 2M* f (§|a| 2 - 1) IS | 



d <n M| (l«| 2 + l)* 3 V Oa| 2 + l) 5 



(6.19) 



Proof. That each solves (|6. 16[) follows by changing variables in (|1.1[) . The initial conditions 
follow from (|5.3[) and (|6.5[) . ■ 

The following lemma will be used in the following to approximate the functions H n {') by 
polynomials P(-;M n ,/3 n ). 

Lemma 6.8 For each n let "H n (Cn) solve \6. 16\) - ![6.19\) . Then, for every e > there exists % 6N 
such that for all n > hq the estimates 

Hn{Cn) - ' (6 - 20) 



| H„ (C«) - P(C« ; M„ , /3„ ) | < e | P(C„ ; M n , p n ) I 



dHn(Cn) dP(C„;M n ,/3 n ) 



d( n d( n 
d 2 U n {Cn) d 2 P(( n ;M n ,p n ) 



dQ 



hold in ^ > and as long as 



dQ 



< £ 



< e 



dP(C„;M n ,j3 n ) 



d(n 

d 2 P(t;M ni /3 n ) 



d( 2 



2Co(M„,/3„) 



> e. 



Proof. The proof of this result is a standard bootstrap argument similar to the ones that has 
been used repeatedly. The idea is that the initial conditions (|6.18[) and (|6.19[) tend to the ones for 
P as n —> oo, sec (|6.11[) . Also the term 5 n /(H n ) 3 in (|6.16[) is negligible if n is sufficiently large 

since S n —> as n —> oo (observe that $(r*) = S^ 1 ). On the other hand, the term Rn((n', M„) can 

i 

be approximated by (1 — Mn( n ) 2 as n — > oo. The resulting limiting equation is then (|6.10j) and 
the values of T-L n and of its derivatives can be approximated at Cn = by those of P(-; M n , p n ) and 
its derivatives there. The difference between H n (Cn) and -P(Cn! M„, p n ) can then be approximated 
arguing as in, for example, Lemma 13. II as well as in Lemma 4.3 of [7J. We observe that, as in [7J, 
(|6.2Up implies upon integration of (|6.16[) . that a condition on Co(M„,/3„) — Cn of the form 



(o(M n , f3 n ) - („ 
2C (M n ,j3 n ) 



loE 



(o(M ni P n ) - ( n 

2Co(M„,/3„) 



l<= 



* 1 2 



must be satisfied for n large enough. Then for every e we can choose hq large enough to obtain 
that the e is larger than the solution of S na \ log(x)| + Co a /IC no | 2 = an( i that the initial data 
are close enough to those of P. ■ 

We now prove that /?„ ~ /3*(Af„) as n — > oo. The idea is to use the fact that the derivative 
of the approximating polynomial P(£ n ; M n , j3 n ) at Co(M„,/3„) is of order one (and negative) by 
Lcmma l6.5l Then, we can use Lemma lC.ll of Appendix [Cl that gives the behaviour in the boundary 
layer where TL„ becomes small (near £ = £™")), to conclude that the next polynomial in the outer 
region is close to one having a double zero in the matching region near £ = C™"). 



Lemma 6.9 Suppose that the $ and its derivatives satisfy the assumptions of Provosition \57l\ and 
let {t*} be the sequence found there. Let the sequences {M n } and {p n } be defined by means 
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of J6'.^[ j and \6.5\) , and let the functions /3*(M n ) be as in i6.9\) and the sequence of functions T-L n 
be given by W.15\) . Then, for any e > 0, there exists a L = L(e) > and a no large enough such 
that if 3>(t*) > L and n > no then 

\Pn-P*{M n )\<e\p*{M n )\ . (6.21) 

Also, for all n > no 

H n (Cn) * -K n (Cn ~ Co{M n )) as Cn^(Co(M n ))- (6.22) 

mth K n = - dP ^J M ^>0 (6.23) 

and 

«n(Cn) ^{(n - Co(M„)) 2 as Cn -> (Co(A/«)) + , (6.24) 

On 

where D n is proportional to by a constant of order one, and there exist ol\ and oli € R 
independent of M n such that < ot\ < if n < CJ2- 

Moreover, for all n > uq there exists Eq > smaZ/ enough an £ cr it,n G (CP*™ — £ Oj£™ m + £ o) 

Co(M„) = -^ T f%^-l) and C.(M„-i) = ^- f %^ - l) (6.25) 

and 

Wn-l(Cn-l) - r„_i(C„-l " C*(^n-l)) 2 M Cn-1 "> (C(M n _i)) + (6.26) 

^ rB _ 1= ^ 2 ^y^) >0 , (6 . 27) 



2 < 



Proof. Suppose that n is very large. We apply Lemma 16.81 for n, thus starting £„ = or at 
£ = £*. It then follows that we can approximate H n by the polynomial P(Cn] M n ,f3 n ) in intervals 
of the form ( n G [0, £o(M n , /?„) — ei] with £i > small but fixed and n large enough. This in 
particular implies (|6.22p and (|6.23[) . The fact that K n is bounded from above and below follows 
from Lcmma [6.5l Using then Lcmma lC.ll we obtain that % n (C) can be approximated as a quadratic 
polynomial for £„ = (o(M n ,/3„) + £i- This implies 



A K 5 

H n (( n )-—^(( n -(o(M n )) 2 



< £ -^(( n -C (M n )f, 

On 



for some £2 > small enough and A of order one, and thus (|6.24p follows. We can then replace 
the variables £„ by Cn-i an d H n by H n -\ using ()6. 15|) and applying again Lemma 16.81 we can 
then approximate the function 'Hn-i(C) by one polynomial which has a double root at the value 
of Cn-i corresponding to £„ = Co{M n ,f3 n ) +£%. This implies (|6.24[) . but also (|6.2ip follows by the 
definition of /?», and therefore also (j6.25p . (|6.26p and (|6.27p follow. ■ 



6.3 Proof of propositions 16.11 and 16.21 

In order to prove the propositions we derive information from Lemma 16.91 in the matching region 
around £™ m for n large enough. Let us then assume that $ satisfies the assumptions of Propo- 
sition [ITTJ so that the sequence (|5.3p is well-defined. Let H(£) be defined by means of (|2.12l) and 
(I2.13p . and satisfies (|l.ip . Let the sequence of functions H n be defined by (|6.15l) . so that, by 
Lemma 16771 each such function satisfies f|6. 16[) with initial conditions f|6. 18[) and (|6.19l) . 

We further assume in the following that the approximating polynomials have j3„ = /3*(M„) 
thus they are as described in Lemma 16.41 and we drop the dependency on (3 n in the notation. 

Lemma 16.91 (|6.25p and the definition of the variables Z n give 

Zo(M n )C = ^(M„_x)C-i • (6.28) 
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Remark 6.10 We can now argue that {M n } is a bounded sequence. Indeed, if M n is very large, 
\B. 1 6\ ) implies that Zo(M n ) is very negative. This would imply, using \6. 28\) that ^-1 > (notice 
that Z*{M n _i) > by definition). However, we cannot have > for large n, because 

C -> -oo. 

Using now the definition of £„ (see (|6.15|l ) and (|6.25|) we can compute (£„ — (o(M n j) and 

(Cn-l - C*i M n-l)) to get 

MSCiCn ~ Co(M„)) = MtlC-l(Cn-l ~ UMn-l)) ■ (6.29) 

On the other hand, the definition of the sequence H n (in (|6.15l) ) gives 

\C\ 5 M n H n ((n) = \C-i\ 5 M n -iH n -i«n-i) ■ (6.30) 

We then change variables according to (|6.29|) and (|6.30[) in (|6.24[) in order to write it in terms of 
the variables Cn-i and W„_i. This gives the asymptotic formula, for n large enough, 

2 

Wn-l(Cn-l) - |3 M ^(C»-l " C*(M„-!)) 2 , Cn-l H- (C*(M«-l)) + ■ 

But comparing this to (|6.26[) implies that the approximation 

Y n ^Ml_x\C n -i? = D n \C\ 20 Mj 

is valid for n large enough (here we have also used (|6.17[) ). 

Using the fact that M n is bounded, and also that D n can be estimated from above and below 
by a constant independent on n, we obtain: 

Ci\C\™MJ < T^M^C-il 3 < C 2 \Cn\ 2 °Ml 
for < d < C 2 . Using (|63]t. then 



l£*l^ 10 l If* I 10 

\C-l\~ l£-ll T 

Using now (|6.28p as well as the fact that Zo(M n ) and Z*(M„_i) are bounded from above 
and below for M n and M„_i bounded (cf. ()B.9j) . (|B. 15[) ) . we obtain, for different C\ and C 2 if 
necessary, that 

Using now f|6 . 13[) for n — 1 we can estimate T„_i from above and below by positive constants 
independent on n. Then: 

$(C-i) * cm:)) 10 (6.31) 

for some C > 0. 

We are now in the position of proving Proposition 16. II 
Proof of Proposition [6TT1 Due to Lemma \6. 91 we can assume that (|6.21[) holds for n large. In 
a similar fashion as in [7J we can make rigorous the argument outlined above by combining the 
lemmas [C . 1 1 and \6 . 81 and prove indeed that (|6.31|) hold. Since by hypothesis $(£„ ) — > 00 as n — > 00 
it then follows that $(£*) < $(£*_i) and this gives Proposition 16. II ■ 

For each n, let now r™ m be the value of t at which <!> reaches the minimum in the interval 
( T n> T n-i) as defined in ()6.2[) . We can now prove Proposition 16.21 

Proof of Proposition [6T2l As before, we only give the formal steps of the proof and refer to [7J 
for details. Let £™ in = £(r™") be defined by means of i[2X2]) . Then, by LemmaEH the fact that 
S n <C 1 and Lemma IC. II we can write, to leading order for n large enough, 
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and that, 

= l^ iB | S |£| 5 M n Wn(Co(Mn)) 

It is clear that for n large enough $(r™ m ) approaches by Lemma IB~51 and employing the scaling 



(that is analogous to the one used in the proof of Lemma IC.1[) , give that the following is a valid 
approximation 



Sri 



1 



(<f(r*)) : 



■/i n (0) 



for n large enough. Here we also use that de definition of M n (see 
(see (|6.17[) ). Thus for n large enough one also has, by (|6.31[) . that 



and that S n = ($(r*)" d ) 



$(r™?) 



C III I II 

Sn-1 



Sn- 



C 



(<f(T*)) 2 ° 



hn-l(O) 



for some order one constant C > 0. 

We finally observe that the values /i„(0) are of order one if n is large enough, by Lemma [C. 11 
Also, we can approximately write the quotients |CT m |/ICil = Z (M n ). But each Z (M n ) is an 
order one constant, since the sequence M n is uniformly bounded. Then, since $(r*) — > oo, we 
have that $(t™") — > 0. Moreover, there exists a constant C > (different from the one above) 
such that, for n large enough, 

$(c™) ~ c(<&(Tr")) 10 , 

and thus (|6.2[) follows. ■ 



7 Convergence to the equilibrium point j9 



In this section wc finish the proof of Theorem 12.31 First we prove the following 



Proposition 7.1 Suppose that (<&(t), W(t), *(r), 0(r)) is a solution of h2.19\) - \2.22}) as found in 
Proposition \4-l\ and satisfying \5.1\) . Then, there exist positive constants C\ and C2 depending 
only on a such that 



lim sup $(t) + 



d$(r) 
dr 



rf 2 $(r) 
dr 2 



<C X , 



(7.1) 



and lim inf $(r)>C 2 >0. (7.2) 

Proof. We recall that the solutions found in Proposition 14 . 1 1 are defined for all t 6 (—00, 00) and 
satisfy lim r _ >00 ($(r), W(r), #(r), 0(t)) = (1, 0, 0, tt/2). Moreover, lim T _ ) ._ 00 0(t) = -§, $(r) > 
for any r G (— 00,00) and lim inf 1-^-00 $(t) < 00. We now claim that (|7.1[) holds for some 
Ci > 0. Indeed, otherwise, due to Proposition 15.11 and Proposition 16.11 there would exist a 
sequence {r*} such that lirrin^oo r* = —00 and lim n _ ! . 0O $(t*) = 00 but such that there exits no 
with $(r*_ x ) > $(t*) for all n > n . Then, since limsup^,^ $(r*) = 00 (cf. (|5.4p) it follows 
that $(t* ) = 00, this yields a contradiction and, therefore, (|7.ip is satisfied. 

Suppose now that (|7.2[) is not satisfied, then in particular this implies, probably taking a 
subsequence, that lim ra _ ! . 00 <I>(T™ n ) = but this contradicts (|6.3|) of Proposition 16.21 ■ 

We can now finish the proof of the main result. 
Proof of Theorem I2T31 Due to Proposition gj] there exists a solution of (j23T)| - (j2~2"2"l) de- 
fined for all r € (-00,00) such that lim T _ ) . 00 ($(T), W(t), #(t), (9(r)) = (1, 0, 0, tt/2). Moreover, 
lim T _ i ._ 00 0(t) = — w, < E , (t) > for any t G (—00,00) and liminf T _»_oo $(r) < 00. Then Proposi- 
tion O gives that {71]) and {72]) hold. 
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We now define a sequence of functions: 

$ n (r) = $(t - n) , W r „(r)=W(T-n), tt n (r) = tt(r - n) n= 1,2,3,... 

Using ()2.19p - (|2.22p . (|7.1j) . (|7.2p . standard compactness arguments and the fact that lim T _>._ o 0{r) = 
— ? we can show that there exists a subsequence {nf\ satisfying lim^oo rij = oo and such that 
(t)j W n - (t), (t))} converges uniformly in compact sets of r to a bounded solution of 
(f2TTj) , say ($00 (t), Woo(t), *oo(r)). Moreover, we have ^(t) > Ci > 0, r e (-00,00). Due 
to Proposition 12.11 it follows that ($oo(r), Woo(t), 4' 00 ('r)) is close to P s if r < and |r| is 
large enough. Using the Stable Manifold Theorem it then follows that ( < & o( T ), W / oo( T ), ^00 (t)) 
is contained in the unstable manifold of P s . However, due to Proposition \2A]f v) it follows 
that the only bounded trajectory contained in the unstable manifold of P s is the is the crit- 
ical point itself, thus (<&oo(t), Woo(t), ^oo(t)) = P s - This implies that the sequence Pj = 
($7ij (0), W nj (0), ^> nj {0)),0 n . (0)) converges to the equilibrium p_ as j — > 00. Therefore the points 
Pj are contained in the centre-unstable manifold of p-, whence lim T ^_ oc (<I>(T), W(t), ^(t), 6(t)) = 
P- and the result follows. ■ 
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Appendix 

A Analysis of the solutions of (11.111) 



We recall here the results concerning the following equation 

d 3 $ 1 

(cf. fHH])) that have been shown in [TJ. For simplicity, we henceforth use the same notation for 
the dependent and independent variables as for (|1.5[) . The following holds. 

Theorem A.l There exists a unique solution of l(A.l\) with the matching condition: 

<&(r) ~ —Kt + o(l) as t -> -00 . 
Moreover, the asymptotics o/$(r) for large r is given by: 

$(t) ~ Ft 2 as r -> 00 for some T > . (A.2) 
Finally, the exists a unique solution of HA . 2\l with matching condition 

$(r) ~ —Kt + o(l) as r — > 00 for some K > . 
It also satisfies that there exists a finite t* such that 

$(t)^0 as t^(t») + . (A.3) 
All other solutions satisfy \A . ty) for increasing t, and, for decreasing t , either \A.cfy or 

$(r) ~ Tt 2 as t — > —00 for some T > 

holds. 

The proof of Theorem IA.1I is done by a series of lemmas. The crucial step is to apply the 
transformation 

d$ _i d 2 $ a . .4, .... 

— = $ 3 u ___ = $ 3 w d T = $ 3 dz ( A.4) 

or (XT 
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that reduces (|A.1I) to the system 



dz 



du 

dz 



dv 
dz 



1 + —u v . 



(A.5) 



The lemmas then give the behaviour of the corresponding trajectories and are given below for 
reference. 

The last two equations in (|A.5[) can be studied independently by means of a phase-plane 
analysis. The isoclines of this system are Ti = {(u,v) : v + -|u 2 = 0}, that has du/dz = 0, and 
r2 = : 1 + f^it = 0}, that has dv/dz = 0. The only critical point is p e = (u e ,v e ) = 

((9/5)3 , —(l/3)(9/5)^), and linearisation gives two complex eigenvalues with positive real part, 
namely A± = (l/2)(l/15)s(7 ± Vni). 

We distinguish five regions, R\ to R5, in the phase plane that are separated by the isoclines. 
These are depicted in Figure [T] where the direction field is also shown. 




Figure 1: Phase portrait associated to (|A.5j) showing the direction field. The thick solid lines 
represent the isoclines and the dashed ones the separatrices. 

Standard arguments imply that any orbit on the phase plane eventually crosses the isoclines 
into the region R4 forwardly in z. We recall that R4 = {(u, v) : —u 2 /'i < v < — 3/(5«) if u < 
, v > max(— u 2 /3, — 3/5u) if u > 0} and the field in it satisfies du/dv > 0. If, however, an orbit 
has (it, v) € i?4 at some value of z, it is possible to discern from which of the regions is coming 
from for smaller values of z by identifying the separatrices of the system. We have the following 
result. 

Lemma A. 2 (Separatrices) (i) There exists a unique orbit v = v(u) in the phase plane as- 
sociated to system \ A . 5\) that is contained in R4 for all u E M. Moreover, v(u) has the 
following asymptotic behaviour 

v(u) = — + O (u^\ as u — > 00 , v(u) = — — (1 + o(l)) as u — > — 00 
2 V / 2u 

(ii) There exists a unique orbit v = v(u) in the phase plane associated to system \A . 5\) that has 
the following asymptotic behaviour 

v(u) = — — (1 + o (1)) as u — > +00 
2u 

and 

\{u,v) - (u e ,«e)l < Ce Re{x+)z as z^-00. (A.6) 
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Lemma A. 3 All orbits associated to solutions of \A . 5\) enter P4. Those that are below the sep- 
aratrix v = v(u) come from the critical point p e as z — > — oo and satisfy t A. 6\) . All other orbits, 
except for v, come from the region R$. 

B Analysis of the polynomial solutions of (11.101) 

Wc identify the functions giving the leading order behaviour of the solutions of if |£| 2 / 3 P(£) 
becomes large on a bounded interval around some negative value £ = £* that gives a large local 
maximum of $ (after the change of variables (|2.12[1 ). 

Due to (|2.13|) . these approximating functions, if £ — > — oo and |£| 2 / 3 Pf is large, are 

ho = \^h(o, (b.i) 

where solves (fTTTUj) . 

Let £ = £* < be such that •!>(£*) is a local maximum. Wc shall derive in the following some 
properties of these functions, such us conditions to ensure that they have a double zero at a point 
£mm ^ £* _ g rs ^. we neec [ t normalise them in an adequate way. 

We can readily integrate (|1.10j) with given initial conditions at £ = £,*■ This gives a family 
of fifth order polynomials that depend on £*, H(£,*) and the first and second derivatives of H 
evaluated at £ = £*. Using that d$>(£*)/d£ = we can write 

thus we can eliminate one parameter and obtain that 

B(o = - ({ ~f )3 (g - o 2 + - n + nn 2 ) 

+ g(r) ( 1 + Sl_n) + ^f)K-r)\ (B . 2) 

We want to characterise the functions l> in terms of the parameters £*, $(£*) and d 2< I>(£*)/<i£ 2 , 
however. In order to do that we first compute 

$(d - iri^r), -^H - -yirr^r) + in 3 d ^ 

and then define a family of polynomials P that normalises H as follows: 

H(0 = \e\ 5 P(Z;M,(3), Z=^, 

thus P solves 

d 3 P 

= < a3 > 

and the parameters M and /3 arc defined by 



1 d 2 $C£*l 

M = P(l;M,/3) and P = „ J . (B.4) 

Then, (|B.1|) implies that 

$(0 = ri ¥ (^) i ^(2;M,^). (B.5) 
Changing to these variables in (|B.2|) we obtain 

P(Z; M,P) = (Z ~ 1)3 (Z 2 + 3Z + 6) + ^(5Z 2 - 16Z + 20) + - l) 2 . (B.6) 
60 9 2 

We point out that Z has the opposite sign as £ and thus we shall be interested in the region 
Z > 1, or that one ahead of the local maximum of (Z) 2 / 3 P '. We next show the existence, for each 
M > 0, of a unique /3*(M) such that P(Z; M,/3) has a unique and double root in {Z > 1}. 
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Lemma B.l Let P(Z; M, 13) be as in HB.6\) with M > as in \B.$ - Then, there exists a /3*(M) < 
such that 



(i) If (3 > /3*(M), P{Z;M,(3) is strictly positive in Z > 1. 

(ii) If f3 < /3*(A/) there exists a unique Z r = Z r (M,(3) > 1 such that P{Z;M,j3) > in 
1 < Z < Z r , P{Z r ; M, (3) = and dP(Z r ; M, (3)/alZ < 0. 

(Hi) If 13 = 0*(M) there exists Z*{M) > 1 that is a double zero ofP{Z; M, (3), and P(Z; A/, (3) > 
for\<Z+Z*(M). 

Proof. The polynomial P(Z; Af, (3) is monotonically increasing in (3. Since (§Z 2 — ifZ + ^p) > 0, 
it follows that P{Z; M : 0) > if Z > 1 and /3 > 0. On the other hand given M > and Z > 1, 
P(Z; M, (3) < for negative values of /3 if |/3| is large enough. 
For each M we define 



/3*(M) = sup {f3 < : P{Z; M, (3) < for some Z > 1} . (B.7) 

Then, by continuity in /3 there exists a Z*(M) such that P(Z*(M); M, /3*(M)) = 0. Moreover, 
dP(Z*{M)\M,P*(M))/dZ = 0, since otherwise there would exists /3 > /3*(M) and Z > 1 such 
that P(Z; M, (3) < 0. 

We now claim that d 2 P(Z*(M); M, f3*(AI))/dZ 2 ^ 0. Indeed, otherwise the definition of /8*(M) 
would imply that Z = Z*(M) is a zero of P(Z; M, /3*(M)) of order four and, in particular, that 
d 3 P(Z*(Af); M, /3*(M))/dZ 3 = 0. But P solves (HU, thus 

> for any (B.8) 

Suppose that j3 < /3*(M). The monotonicity of P(Z; M, /3) in /3, combined with the fact that 
P(1;M, (3) = M > and limz^oo P(Z; A/, /3) = oo imply that there exists at least two zeros of 
P(Z;M,f3) in {Z > 1}. Actually, there are exactly two zeros of P(Z; A/, (3) in {Z > 1}. For 
otherwise, there would be four zeros in {Z > 1} counting multiplicities, and this would imply, by 
Rolle's Theorem, the existence of three zeros of dP(Z; M, (3)/dZ in {Z > 1}, and iterating the 
argument, also the existence of two zeros of d 2 P(Z; M, (3)/dZ 2 in {Z > 1} and at least one zero 
of d 3 P(Z;M,/3)/dZ 3 in {Z > 1}. But this contradicts (fBTSj) . 

Therefore, P(Z; Af, 0) has exactly two zeros in {Z > 1}. The smallest of which, Z r , satisfies 
dP(Z r ;M, 13) /dZ < (by continuity since M > 0). ■ 

We now compute the asymptotic behaviour of /3*(M) in the limits M — > and M — > oo. 

Lemma B.2 Por every M > let /3*(Af) and Z*(M) 6e as m Lemma \B.l\ Then, they satisfy 
that 1 < Z*(M) < 4 and /3*(Af) < 0, and /iave i/ie following asymptotic behaviour: 



Z*(Af)~ l + (12Af)3 , /3,(M) 
d 2 P(Z*;M,(3*(M)) /3Af\^ 

dZ 2 



(B.9) 



as A/ -> 



Q 71 T 

Z,(M)->4, &(M) -, 

d 2 P(Z*; A/, f3*(M)) 2M 

5z3 r as 



(B.10) 
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Proof. In order to determine Z*(M) and 8*(M) we have to solve the system that results from 
imposing that P(Z; M, B) has a double zero: 

(Z ~ Q 1)3 (Z 2 + 3Z + 6) + y (5Z 2 - 16Z + 20) + |(Z - l) 2 = , (B.ll) 

V 10 ' {Z 2 + 2Z + 3) + — (5Z - 8) + B(Z - 1) = . (B.12) 
12 9 

Subtracting the second equation of (|B.12j) multiplied by (Z — l)/2 to (jB.llj) we obtain 

(Z - 1) 3 (3Z 2 + 4Z + 3) = 40(4 - Z)M , (B.13) 

the solution of which gives Z*(M), the position of the double zero. Now B*(M) can be computed 
from either (|B~TT|) or (|BTl2|) by substituting Z = Z*(M). It is clear that Z*(M) G (1,4) for any 
M > 0, since the left-hand side of (|B.13|) is positive for Z > 1. 

It then follows from (|BTT3|) that, if M -> 0, (Z* (M) - 1) is of order Ma and that, using (|B.lip . 
B*(M) behaves like Ms. On the other hand, if M — > oo, we obtain that Z*(M) — > 4~ and that 
/3*(M) is of order M. The precise asymptotic behaviours stated in (|B.9|) and in (|B.10|) follow easily 
from (|B.13|) and (jB.llj) by using the leading order behaviour of Z*(M) in both limits M — > and 
M -> oo. ■ 

Finally, we identify the largest root of P in the region Z < 1: 
Lemma B.3 For all M > and /3 £ K i/ie following value is well defined 

Z (M, B) := max{Z < 1 : P(Z; M, /3) = 0} . (B.14) 
If B = 8*{M) then, setting Z (M) := Z (M, B*(M)) 

Z 0( M)~1-^, dP( !° ;M) ~m § Mf a. M^0, (B.15) 



dZ V 2 , 
and 

„ /20M\* dP(Z ;M) l/20\',,4 , r „, 

Z (M)~-(— J , 3( T J a, M-oo. (B.16) 

Moreover, in this case, Zq(M) is the only real root in Z < 1. 



The proof follows by continuity, Lemma IB.ll and the Implicit Function Theorem. 
We need to derive some information concerning the derivative of the polynomial P(Z; M, 8) 
at Z = Z (M,/3). 

Lemma B.4 The value Z (M,f3) is the unique root of P(Z; M, 8) in {Z < 1}. Moreover, if 
Zq(M,/3) > — 1, there exists a positive cq independent of M and f3 such that 

dP(Z (M,B);M,B) r 2 , rl 

i-^ — '—^ — '-t-L > c max{M3,M}. B.17) 

dZ 

Proof. Using ([B3|) . and differentiating P(Z;M,B) then 

fM|^)4 (Z 3_ 1)+(T+ , 

and we obtain that P(Z; M, 8) is strictly concave for Z < 1 if 10M/9 + /3 < 0. On the other hand, 
if WM/9 + B > we obtain that P{Z;M,B) is strictly concave for Z < Z := -(l-3(10M/9+/3))s 
and strictly convex if Z > Z. Then, since dP(l; M, B)/dZ = — 2M/3 it follows that there exists a 
Z(M, /?) < 1 such that dP(Z; M, 8) jdZ > if Z < Z(M, /3) and dP(Z; M, /3) /dZ < if Z(M, /?) < 
Z < 1. Using now that P(1;M, B) = M > it then follows that P(Z;M,B) attains a positive 
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maximum in the interval Z E (—00, 1) at Z = Z(M,/3). Therefore, since lim^^-oo P(Z; M,/3) = 
—00 and the concavity, the function P(-; M, /3) has a unique zero Zq(M, j3) £ (—00, 1). 

It is easy to prove that for M bounded dP{Z]M,j3)/dZ > Kq > uniformly for /3 < 0, 
uniformly on Z £ (—1, 1) (differentiating (|B.6[) V In fact, one can actually show that for M small 
enough there exists a positive constant Co such that dP(Z; M, f3)/dZ > cqAIs if Z £ (—1,1). 

In order to obtain (|B.17[) for large M we first consider the case \(3\ < EqM with £0 sufficiently 
small, and then (|B.17[) follows easily Z £ (—1,1) Suppose then that (3 < —sqM. Then, if M is 
large we can approximate P(Z; M, f3) and its derivatives in the interval Z £ (— 1, 1) by 



M 



g -(5Z 2 - 16Z + 20) + |(Z - l) 2 = MW M (Z) 

and its derivatives respectively. Observe that then Wm{Z) is a quadratic polynomial with bounded 
coefficients satisfying Wm(1) = 1. Suppose that Z = Z (M,(3) e (—1, 1). Notice that we cannot 
have any other zero of Wm(Z) in the region Z < 1 because then, by continuity, P(Z; M, j3) would 
have more than one zero in the domain {Z < 1} and this would contradict the statement above. 
Then, using also that Wm(1) = 1, we can write 

{Z-Z a ){Z-Z x ) 



W M {Z) 



z, 



with (x-z^i-Zx) bounded and Z\ > 1. Therefore min{(l 



0--z o )(i 

Z ), (Z l 

for large M. It then follows that (Z\ - Z ) > 2K\. Then, since Zq > — 1 

dW M (Z ) _ (Zo-Zj (Z 1 -Z ) ^1(Z X -Z 

dZ 

And the result follows. 



> 



(1 - Z )(l - Zi) (1 - Z„)(Zi - 1) " 2 (Zi - 1) 



1)} > i^i > 0, uniformly 



> ci > . 




Figure 2: Schematic depiction of the polynomials P(Z;M,j3) for fixed M and different values of 
/?. The solid line represents a polynomial with /3 = /3*, the dashed line one with /3 > /3» and the 
dashed-dotted line one with (3 < [3*. The figure also reflects the fact that dP(l; M, j3)/dZ < for 
any j3 < and M > 0. 



C Analysis of the bouncing region 



We seek to reformulate the results concerning the region of very small H (or $) that can be 
approximated by p. lip , that were obtained in [7]. The result here is more general, namely the 
class of equations under consideration is 

■ R(Q = ° 5 > (C.l) 



dC 3 



where, if C belongs to some given bounded interval, then 



\R(0\ 



-(C) 
d( v ^ 



<c . 



(C.2) 
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Lemma C.l Suppose that the equation \C. 1\) is satisfied in an interval C € (Co — £i s Co + £1) for 
some ei > independent of S and where R satisfies HC.ty) . Suppose that: 



\n(o+m-co)\<s2K\c-c \ 



<M 

d( 



(0+K 



< e 2 K 



d 2 H 



d( 2 



(0 



< £2 



(C.3) 



for C € (Co — £i,Co — £ i/2), K > c > cmrf some e 2 > 0. For any e 3 > 0, there exists e > 
independent of 5 and K , but in general depending on cq, and a real constant A > such that if 
£ i + £2 < £o i/iere exists a 5q = 6o(si, £3, cq) > such that for 5 < 5q we have: 



no 

<£H. 



(0 



AK 5 
6 

2AK 5 



d 2 H 



dC, 2 



(0 



(C-C0) 2 
(C-Co) 

2AK 5 



< 



< 



5 

£ 3 K 5 
S 

£ 3 K 5 



Co) 2 - 

(C-Co), 



(C.4) 



/orCe (Co+£i/2,Co+£i). 



Proof. This result can be adapted from the results proved in [7]. More precisely, the result follows 
from arguing as in the proof of the lemmas 4.4, 4.5, 4.6 and 4.7 of this article. We recall the main 
ideas of the argument here. 

We introduce the variables h(s) = K 3 H(Q/5 and s = K 4 (( — Co)/^ and obtain that h satisfies 
the equation 

f 5 2 \ 1 

5* — 



cPh 
ds^ 



h 3 



(C.5) 



with the matching condition 



h(s) 



We can then reformulate (|C.5|) by the transformation (|A.4[) of Appendix [A] thus, to leading order 
of approximation and due to the assumption (|C.3p , the solution of the resulting perturbation of 
(IA.5I) follows the scparatrix v(u) (see Lemma |A.2[) . A key point in the argument is that the system 
of ODEs (|A.5|) is integrated forward, the direction for which the separatrix is stable. This allows 
As 2 as s 



to prove that h(s) ~ As 2 as s — > 00 for some A > (see Theorem IA.1[) . thus the asymptotics of 
the resulting solution for C > Co can be described by means of (|C.3p . ■ 
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